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ABSTRACT 

(Distribution  Limitation  Statement  B) 


Coupling  of  transients  to  aboveground  transmission  lines  is  analyzed  and 
measured  experimentally.  Measurements  were  made  using  scale  models  approxi¬ 
mately  1-foot  high  to  represent  power-  transmission  lines.  The  coupling  to  the 
line  and  the  propagation  characteristics  of  the  line  were  measured.  The 
effects  of  bends  and  junctions  on  the  transmitted  signals  were  measured  and 
analyzed.  Periodically,  grounded  lines  and  mult"' conductor  lines  were  also 
examined.  The  variation  of  the  induced- voltage  waveforms  with  soil  conductivity, 
line  height,  and  pulse  width  were  calculated. 
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SECTION  I 


INraODUCTION 


A.  Background 

1.  Previous  Work 

The  coupling  of  tii®  E5iP  wave  to  power  transaission  lines  was 
initially  analyzed  using  a  low-frequency  version  of  Sunde's  theory,  ~ 
which  was  valid  for  a  long,  straight  transaission  line  whose  height  was 
scalier  than  the  longest  wavelength  of  interest.  The  analysis  of  trans- 
aission  lines  using  this  approach  was  reported  in  Technical  Meaorandua 
22  for  this  contract."^  Because  of  its  apparent  lisitation  to  the  low 
frequency  portion  of  the  EJI?  spcctrua.  however,  two  prograss  were  initi¬ 
ated  to  isprove  the  iindcrstaralii:^  of  the  coupling  to  the  transaission 
lines.  The  first  was  a  core  rigorous  analysis  of  the  coupling  to  the 
trar^ission  line  that  was  not  liailed  to  a  height- to-wavelength  ratio 
less  than  one  (hx^r.  <  1) .  The  details  of  this  effort  will  be  reported  in 
a  Separate  doraaent;  however,  sosc  of  tne  ieportant  results  are  contained 
in  .App*--r.3ix  A.  The  second  was  an  experiacntal  dctcralnatlos  of  the  actual 
coupling  to  a  transsission  line  when  the  hcight-to-wavelecgth  ratio  b/\ 
is  not  scalier  than  one.  To  study  coupling  to  a  power  transsisslon  line 
in  such  3  sanner  that  directional  effects,  propagation  characteristics, 
radiation,  and  scattering  effects  could  be  adequately  assessed,  however, 
w-tald  r.>.-qulre  an  extensive  effort,  a  large  test  site,  and  a  large  aobile 
tranSHitter.  The  cost  of  a  thorough  evaluation  of  the  coupling  w^ld  have 
prohibitive,  and  the  results  that  could  have  been  obtained  free  a 
sraple  t*ist  would  have  answered  fewer  questions  than  they  left  unanswered. 
Because  of  the  licltations  in  the  practical  ability  to  experisentally 
exanine  the  coupling  to  full  scale  power  transaission  lines,  at^  bccauss 
it  was  anticipated  that  the  results  of  the  rigorois  analysis  of  the 
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coupling  Bight  not  be  available  soon  enough  to  meet  the  icKOdiate  require¬ 
ments  of  the  Air  Force,  it  was  propose^;  to  use  scale  sodellis^  of  the 
power  lines  to  study  fhe  coupling  expcri^-ntally. 

2.  •  Use  of  Scale  Models 

The  use  of  scale  aodels  to  obtain  ei^neerii^  solutions  to 
coaplicated  electroeagnetic  probleas  is  an  old  but  powerful  techniip^. 

It  is,  in  essence,  an  analog  coaputer  technique  used  for  probleas  that 
arc  very  difficult  to  solve  accurately  Cor  with  confidence)  by  com-ec- 
t tonal  approaches  used  in  oatheoatics  and  physics  either  because  they 
entail  coaplicated  geoaetries  or  because  the  relative  ioportance  of  the 
various  factors  thought  to  influence  the  solution  is  unhnoMm.  These 
probleas  can  be  solved,  however,  by  constructing  electrooagnetic  oodels 
of  the  structures,  surroundings,  and  excitation  sources,  and  seasuring 
the  solution.  A  classic  example  of  the  use  of  scale  oodels  is  the  use 
of  electrolytic  tanks  to  nap  the  fields  or  equipotential  surfaces  about 
irregular-shaped  bodies  that  are  (or  were)  difficult  to  analyze.  A  oorc 
current  cxaoplc  is  the  use  of  models  to  study  aircraft  antcniai  perfor- 
oancc.  Since  the  typical  airfrase  is  such  too  cfxapllcated  to  persit  an 
analytical  dctcraination  of  antenna  patterns  and  icpcdanccs  with  suffi¬ 
cient  cconosy  (even  with  aodem  computation  sethods)  to  justify  aircraft 
antenna  design  based  on  analytical  techniques  alone,  aircraft  (and  other 
vchioilar)  antenna  design  has,  for  the  last  two  decades,  relied  on  the 
scale  Bodcl  sethod  of  dctcrainis^  radiation  patterr>s. 

Elcctrosagnctic  scale  iKKiclling  is  useful  prinarlly  because 
the  technique 

•  Is  cconoai<»l. 

•  Allows  zlBCly  detetoination  of  engineering  solutions. 

•  Provides  solutions  fo  probleas  too  difficult  to  solve 
analytically. 


•  Permits  greater  control  to  be  exercised"  over  the  experiment 
under  laboratory  conditions. 

•  Allows  flexibility  to  modify  the  model  to  study  the  effects 
of  changes  in  geometry,  an  advantage  that  is  particularly 
useful  when  the  exact  geometry  is  unknown. 

The  technique  is  particularly  well  suited  to  the  solution  of  problems 
that  are  too  complicated,  because  of  the  geometry  or  parameters  involved, 
to  solve  in  a-  timely  manner  by  analytical  methods,  yet  simple  enough 
structurally  that  a  scale  model  can  be  constructed  quickly  and  economi¬ 
cally.  'Power  transmission  lilies  and  airframes  are  good  examples  of 
structures  that  can  be  constructed  quickly  and  economically.  The  interior 
of  the  aircraft  or  the  structure  inside  a  power  transformer  are  much  more 
difficult  to  build  into  scale  models;  hence  they  are  not  good  candidates 
for  analysis  using  scale  models  (at  leas,t  not  with  the  same  models  used 
to  study  the  power  line  and  the  airframe),. 

3.  Experimental  Verification  of  Coupling  to  Transmission  Line 

The  scale  model  facility  used  to  study  the  common  mode  coupling 
to  the  transmission  line  consisted  of  an  area  of  soil  that  was  salted  to 
increase  its  conductivity  to  about  0.2  mho/m.  The  test  area  was  illumi¬ 
nated  with  a  pulse  radiated  by  a  dipole  antenna.  The  radiated  pulse  was 
essentially  a  two-exponential  pulse  with  a  rise  time  of  less  than  0.3  ns, 
and  a  decay  time  that  could  be  varied  from  a  fraction  of  a  nanosecond  to 
about  35  ns,  so  that  responses  ranging  from  the  impulse  response  to  the 
step-function  response  could  be  obtained.  Both  vertical  .polarization 
and  horizontal  polarization  of  the  wave  were  used,  and  azimuth  angles  of 
incidence  from  0  to  180  degrees  relative  to  the  direction  of  the  trans¬ 
mission  line  were  examined.  Elevation  angles  of  0  and  about  30  degrees 
with  respect  to  the  horizontal  were  used  for  most  of  the  experiments. 

The  experiments  were  based  on  a  nominal  scale  of  20:1,  so  that  a  20-foot 
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high  transmission  line  over  soil  of  conductivity  10  mho/m  would  be 
simulated  by  a  1-foot-high  transmission  line  over  salted  earth  with  con¬ 
ductivity  0,2  mho/m.  This  basic  model  could  also  be  interpreted  as  a 
40:1  scale  model  of  a  40-foot- high-  transmission  line  over  soil  of  con- 

_3 

ductivity  5  X  10  mho/m,  however.  In  fact,  in  the  course  of  the  para¬ 
metric  study  of  coupling^  the  model  line  height  was  varied  from  0  to  2 
feet,,  so  that  a  range  of  line  heights  (or  by  .a  different  interpretation, 
a  range  of  soil  conductivities)  was  studied.  The  scaling  considerations 
for  electromagnetic  modeling  are  discussed  in  detail  in  Appendix  B. 

The  results  of  the  scale  model  studies  of  common  mode  coupling 
to  the  transmission  line  indicated  that,  for  the  full  scale  -rise  times 
of  5  to  10  ns,  the  low-frequency  analysis  developed  in  Technical  Report 
b  for  this  contract®  is  adequate  to  describe  all  significant  features  of 
the  transmission  line  response.  This  result  was  somewhat  surprising  in 
view  of  how  far  beyond  the  limits  of  validity  of  the  low-frequency  anal¬ 
ysis  the  spectrum  of  a  pulse  with  5  ns  rise  time  extends.  The  low  fre¬ 
quency  analysis  represents  the  induced  current  so  accurately  because 

(1)  the  field  strength  of  the  wave  scattered  from  the  wire  and  reflected 
back  to  the  wire  by  the  ground  is  small,  so  that  the  current  it  induces 
in  the  wire  is  small  compared  to  that  induced  by  the  direct  wave,  and 

(2)  over  the  pulse  spectrum,  the  "scattering  impedance"  (i.e.,  the  ratio 
of  the  axial  component  of  the  incident  field  to  the  current  induced  in 
the  wire)  is  nearly  equal  to  twice  the  characteristic  impedance  of  the 
wire  over  ground,  so  that  the  induced  current  is  almost  the  same  whether 
the  wire  is  viewed;  in  early  times  as  a  transmission  line  (low-frequency 
model)  or  as  a  scattering  cylinder  in  space.-  (For  much  faster  rise  times 
there  should  be  a  noticeable  difference,  however,  because,  for  frequencies 
considerably  higher  than  those  contained  in  the  piilse  spectrum,  the 
"scattering  impedance"  is  significantly  smaller  than  2Z^,  so  that  the 
induced  current  at  very  early  times  would  be  significantly  larger  than 

that  predicted  by  the  low-frequency  transmission-line  theory.) 
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The  experimental  studies  of  coupling  to  long,  straight  trans¬ 
mission  lines  thus  demonstrated  that  the  hypothesized  coupling  mechanisms 
were  correct  and  ihat  the  effect  of  line  height,  soil  conductivity,  and 
angle  of  Incidence  predicted  by  the  low-frequency  transmission-line  theory 
were  adequate  throughout  the  pulse  spectrum  of  interest  (i.e.,  10  kHz 
to  100  MHz) .  In  addition,  the  effect  of  coupling  to  a  vertical  element 
such  as  a  transformer  ground  wire  or  a  service  entrance  conduit  were 
studied  using  the  scale  model,  and  the  transmission-line  theory  was 
adapted  to  the  analysis  of  coupling  to  these  elements. 

B.  Other  Aspects  of  Coupling  to  Power  Lines 

A  typical  power-distribution  transmission  line  is  not,  however,  a 
long  straight  line.  The  typical  distribution  line  contains  bends  neces¬ 
sary  to  follow  the  right  of  way  and  spur  lines  that  branch  off  from  the 
main  line  to  serve  other  loads  within  the  distribution  area.  In  addition, 
the  transmission  line  contains  three  or  four  conductors  instead  of  the 
single  conductor  assumed  in  the  common-mode  model,  and  in  the  Tour-wire 
wye  system,  the  neutral  wire  is  sometimes- grounded  periodically  along 
the  transmission  Tine.  Since  these  features  of  practical  distribution 
lines  may  affect  the  coupling  to  the  lines  or  the  propagation  of  the 
coupled  signal  to  the  customer's  terminals,  additional  model  studies, 
accompanied  by  analyses,  were  made  to  evaluate  the  effects  of  these 
phenomena. 

In  addition,  the  analysis  of  the  coupling  presented  in  Technical 
Report  5^  was  in  the  frequency  domain,  and  much  of  the  scale  model  work 
was  performed  with  pulses  approximating  a  step  function  because  the  step 
function  response  is  readily  analyzed  and  interpreted.  The  exponen¬ 
tial  pulse  response  is  more  representative  of  the  however.  Now 

that  the  coupling  theory  has  been  validated  using  the  model  studies,  the 


response  of  the  transmission  line  to  exponential  pulses  for  various  line 
heights  and  soil  conductivities  can  be  calculated  to  assess  the  effects 
of  these  parameters  on  the  waveform  and  peak  values  of  the  induced  cur¬ 
rent  or  voltage. 

C.  Purpose  of  this  Report 

The  purpose  of  this  report  is  to  examine  the  effect  of  the  coupling 
parameters  on  the  induced  waveform  and  to  examine  the  effects  of  the 
non-ideal  aspects  of  practical  transmission  lines  on  the  coupling  and 
propagation  of  the  induced  signals.  In  addition^  the  status  of  the  power¬ 
line  coupling  was  to  be  reviewed  to  determine  the  requirements  for  full- 
scale  testing  on  a  power  transmi.ssion  line,  if  such  testing  is  deemed 
necessary. 

D.  Summary  of  Results 

I.  Common-Mode  Coupling  to  the  Transmission  Line 

The  maximum  commoiiTmode  coupling  to  the  transmission  line,  as 
measured  by  the  open-circuit  voltage  or  short-circuit  current  at  the  end 
of  a  long^  straight  line,  is  obtained  for  a  vertically  polarized  wave 
incident  along  the  axis  of  the  transmission  line  and  propagating  toward 
the  terminalSi  For  a  truly  semi-infinite  line  length,  the  peak  voltage 
is  limited  only  by  ihe  attenuation  of  the  transmission  line,  but  for  the 
more  likely  case  of  a  finite  length  of  straight  line  a  few  kilometers 
long,  the  peak  voltage  is  limited  by  line  length.  The  waveform  of  the 
short  circuit  current  induced  at  the  end  of  a  finite  length  of  straight 
line  by  the  vertically  polarized  exponential  wave  incident  along  the  axis 
is  a  short  spike  of  current  having  a  rise  time  approximately  equal  to  the 
rise  time  of  the  incident  field  and  a  fall  time  roughly  one-tenth  that 
of  the  incident  field  (assuming  the  rise  time  of  the  Incident  field  is 
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much  less  than  its  fall  time).  For  other  angles  of  incidence,  the  rise 
time  of  the  open-circuit  voltage  (or  short-circuit  current)  is  a  few 
tenths  of  the  fall  time  of  the  incident  wave,  and  the  fall  time  is  very 
slow  (several  fail  times  of  the  incident  wave) .  In  either  case,  the 
peak  current  or  voltage  is  roughly  inversely  proportional  to  the  square 
root  of  the  soil  conductivity  and  only  weakly  dependent  on  line  height 
(except  when  the  line  is  oyer  very  highly  conducting  soils) . 

-2 

For  average  soil  conductivity  (10  mho/m),  a  lO-meter  line 
height,  andi  a  1-V/m  exponential  pulse  with  a  l-ps  decay  time  incident  on 
a  1000-meter  long  line,  the  peak  induced  current  will  be  about  0.93  A  for 
a  plane  wave  and;  abou-  0.65  A  for  a  spherical  wave  (with  1  V/m  at  the 
end'  closest  to  the  source) .  For  the  lossless  semi-infinite  line  with 
the  plane  1  V/m  wave  incident  at  an  elevation  angle  of  10  degrees,  the 
short-circuit  current  would  be  2.8  A.  The  corresponding  open-circuit 
voltages  for  a  300-ohm  transmission  lino  are  280,  195,  and  840  volts  per 
V/m,  respectively.  It  should  bo  emphasized  that  those  responses  arc  for 
vertically  polarized  waves  incident  at  gr.azing  angles  along  the  axes  of 
long  straight  transmission  lineSi  Whether  this  combination  of  source 
characteristics  and  transmission  line  characteristics  actually  occurs 
frequently  in  practical  situations  of  concern  should  be  considered  before 
such  coupling  characteristics  are  used  as  system  excitation  criteria. 

2.  Propagation  Characteristics  of  Transmission  Line 

The  propagation  characteristics  of  the  common  mode  signal  were 
studied  experimentally  by  observing  the  change  in  the  leading  edge  of  a 
fast-rising  pulse  as  it  propagated  along  the  line,  and  by  measuring  the 
position  and  width  of  the  null  in  vhe  CW  standing  wave  on  the  line  when 
it  was  short-circuited  at  one  end.  in  both  cases,,  the  results  were  found 
to  agree  very  well  with  those  given  by  Sunde^  and  Technical  Report  5.^ 
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It  was  concluded^  therefore,  that  the  dominant  parameters  affecting  the 
propagation  characteristics  of  the  power  transmission  line  were  the  con¬ 
ductor  height  and  the  soil  conductivity  and  that  these  factors  are  ade¬ 
quately  accounted  for  (for  engineering  purposes)  by  the  work  of  Sunde.^ 

3.  Effects  of  Bends  and  junctions  on  Propagation 

A  bend'  in  the  transmission  line  was  found  to  cause  a  small 
reflection  of  the  leading  edge  of  the  incident  current.  The  properties 
of  the  tr'ansmission  line  can  be  represented  by  a  distributed  decrease  in 
the  inductance  and  capacitance  per  unit  length  in  the  vicinity  of  the 
bend,  which  causes  a  local  decrease  in  the  effective  characteristic  imped¬ 
ance  of  the  line.  The  peak  reflected  current  for  sharp  bends  of  90  de¬ 
grees  or  less  is  6  percent  or  less  of  the  incident  current  for  full  scale 
current  rise  times  of  the  order  of  10  ns.  Because  a  little  of  the  lead¬ 
ing  edge  of  the  incident  current  pulse  is  reflected,  the  rise  time  of 
the  current  transmitted  beyond  the  bend  is  slightly  longer,  but  these 
effects  are  quite  small  and  can  be  neglected  without  serious  error  for 
bends  smaller  than  90  degrees.  If  the  radius  of  the  bend  is  increased, 
the  effect  of  the  bend  is  reduced,  so  that  for  a  bend  radius  of  20  feet 
in  a  20-foot-high  line,  the  reflection  is  reduced  from  6  percent  to  2 
percent  for  a  90-dogree  bend. 

The  effect  of  a  branch  in  the  line  formed  by  the  junction  of 
two  load  lines  with  a  source  line  is,  to  a  first  approximation,  repre¬ 
sented  well  by  classical  transmission  line  theory.  When  a  source  line 

of  characteristic  impedance  Z  is  loaded  by  two  identical  transmission 

o 

lines  in  parallel,  the  current  transmitted  to  each  of  the  load  lines  is 
two  thirds  of  the  current  incident  on  the  junction  from  the  source  line, 
and  the  current  reflected  back  down  the  source  line  from  the  junction  is 
one  third  the  Incident  current.  When  the  angle  between  the  two  load  lines 
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is  small,  there  is  some  deviation  fi'om  this  ideal  line  behavior  because 
of  mutual  coupling  between  the  load  linos.  The  current  transmitted  to 
the  load  lines  decreases  when  mutual  coupling. is  significant,  however, 
so  that  the  maximum  transmitted  current  occurs  when  the  lines  are  behaving 
as  ideal  transmission  lines. 

4.  Multiconductor  Transmission  Lines 

The  studies  of  multiconductor  transmission  lines  indicate  that 
differential  currents  whose  magnitude  is  related  to  the  wire  spacing  and 
angle  of  incidence  of  the  wave  ai*e  induced  in  the  transmission  lines. 

The  results  are  generally  consistent  with  those  predicted  by  Harrison 
et  al.,  but  their  work  was  for  two-wire  transmission  lines  in  free  space, 
while  the  cases  of  interest  here  are  for  two  or  more  wires  above  a -ground 
plane.  Because  the  ground-reflected  wave  incoin)letely  cancels  the  direct 
wave,  the  late-time  differential  current  does  not  go  to  zero  as  might  be 
expected  for  two  wires  over  a  perfect  ground.  Instead,  a  residual  dif-s 
ferential  current  of  about  20  percent  of  the  peak  value  continues  to  flow 
at  late  times.  For  typical  line  heights  and  wire  spacings  the  peak  dif¬ 
ferential  current  for  a  45-degree  angle  of  incidence  can  be  of  the  order 
of  one-tenth  of  the  common  mode  current  with  a  perfect  ground  for  hori¬ 
zontal  polarization  incident  at  broadside.  At  smaller  angles  of  incidence 
the  fraction  can  be  larger,  because  the  common  mode  coupling  decreases 
and  the  differential  mode  coupling  increases  as  the  elevation  angle  de¬ 
creases,  The  differential  mode  coupling  tends  to  be  sisuilest  "'hen  the 
common  mode  coupling  is  largest,  however,  so  that  the  maximu;;  .fferential 
cui'rent  is  quite  small  compared  to  the  maximum  common  mode  current.  Thus, 
for  example,  the  ma.ximum  differential  current  lor  1  V/m  incident  on  a 
two-wire  line  with  one-meter  spacing  is  of  the  order  of  6  to  10  mA, 
whereas  maximum  common  mode  currents  of  the  order  of  I  to  3  A  (at  a 
different  angle  of  incidence)  wore  calculated. 


For  a  four-wire  system  with  the  neutral  periodically  grounded, 
however,  the  differential  current  between  the  three  phase  "on^uctors  and 
the  neutral  can  be  almost  as  large  as  the  common  mode  cun  e  '  J,  apparently 
because  the  neutral  is  maintained  at  approximately  grot.i  ;  t 
while  the  phase  conductors  assume  '.pproximately  the  potei  t,  they  would 
have  if  the  neutral  conductor  were  not  present. 

5.  Periodically  Grounded  Conductor 

The  common  mode  propagation  charactex’istics  of  a  periodically 
grounded  conductor  were  found  to  behave  as  a  transmission  line  with  many 
spurs  in  early  times  (before  the  reflections  from  the  bottoms  of  the 
ground  wires  returned  to  the  horizontal  conductor)  and  as  a  conductor  in 
contact  with  the  soil  at  late  times.  In  early  times,  the  transmitted 
current  was  roughly  (2/3)”  times  the  incident  current,  where  n  is  the 
number  of  ground  leads  along  the  line.  At  late  times,  the  series  imped¬ 
ance  per  unit  length  can  be  represented  by  2  /?-.  where  Z  is  the  impedance 

&  S 

of  the  ground  rods  at  the  base  of  the  ground  leads,  and  //  is  the  spacing 
between  the  ground  leads.  At  an  intermediate  frequency  associated  with 
the  line  height  and  the  spacing  between  grounds,  there  was  a  minimum  in 
the  line  attenuation,  and  for  long  sections  of  periodically  grounded  line, 
this  frequency  tended  to  dominate  the  transmitted  current. 


SECTION'  II 


COMON-MODE  CURRENT  AND  VOLTAGlS  INDUCED  IN  TRANSMISSION  LINES 

A.  General 

The  analysis  of  the  common-mode  coupling  to  transmission  lines  is 
described  in  consid'-i’able  detail  in  Technical  Memorandum  22.^  There  the 
detailed  analysis  of  the  coupling  and  experimental  verification  of  the 
analytical  results  are  presented.  In  this  section  the  analytical  pro¬ 
cedures  developed  in  Technical  Memoa*andum  22^  will  be  used  to  calculate 
the  common-mode  voltage  and  current  wavefo]ms  for  various  angles  of  inci¬ 
dence,  soil  conductivities,  line  heights,  and  incident  pulse  durations. 

In  Sections  B  and  C,  the  effect  of  these  parameters  is  investigated 
for  the  semi-infinite  horizontal  conductor.  There  It  is  shown  that  the 
strongest  coupling  occurs  for  vertical  polarization  incident  on  a  poorly 
conducting  soil  at  small  angles  of  incide*’ce.  Because  the  length  of  line 
effective  in  producing  a  current  or  voltsge  waveform  of  duration  t  Is  of 
the  order  of  ct/(l  -  cos  cos  <p),  the  semi-infinite  transmission  line 
model  is  valid  for  small  angles  and  cp  only  if  the  line  is  very  long, 
and  for  \Jt  =  9  =  0  the  line  length  must  be  truly  infinite. 

A  more  reall'-'^ic  model  for  the  transmission  line  when  ^  =  cp  =  0  is 
a  line  of  finite  length  with  the  incident  field  sti'ength  varying  as  1/z 
along  its  length.  This  case  is  examined  in  Section  II-D,  where  it  is  known 
that  the  major  coupling  is  due  to  the  horizontal  electric  field  induced 
in  the  soil.  The  waveform  of  the  induced  curi’ent  or  voltage  is  tlien  a 
function  of  the  incident  pulse  shape  and  the  electrical  pi'operties  of 
the  soil.  The  waveforms  for  this  end-fire  coupling  =  9-0)  are  cal¬ 
culated  for  a  two- exponential  incident  pulse  waveform  for  soil  conduc- 
-1  -3 

tivities  between  10  and  10  mho/m.  As  in  the  case  of  the  semi-infinite 
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line,  the  induced  current  increases  as  the  soil  conductivity  decreases 
(I  ^  constant),  but  whereas  the  current  induced  in  a  low-loss,  semi¬ 
infinite  transmission  line  approaches  infinity  as  the  line  attenuation 
approaches  zero,  the  current  induced  in  a  finite  length  transmission  line 
is  always  finite  even  if  there  is  no  line  attenuation.  The  waveform  of 
the  induced  current  for  the  end-fire  condition  is  also  quite  different 
from  that  for  oblique  incidence  because,  in  the  end-fire  case,  all  of  the 
induced  effects  arrive  at  the  terminals  in  synchronism,  while,  in  the 
case  of  oblique  incidence,  there  is  some  dispersion  of  the  induced  cur^ 
rent. 

The  current  induced  in  a  vertical  element,  such  as  a  ground  wire  or 

service  entrance  conduit,  is  calculated  for  an  exponential  incident  pulse 

in  Section  II-E.  In  Technical  Memorandum  22,®  this  current  was  shown  to  be 

significant  in  comparison  with  the  current  induced  in  a  horizontal  con- 

-2 

ductor  over  a  good  ground  plane.  If  the  soil  conductivity  is  10  mho/a 
or  less,  however,  the  current  induced  in  the  riser  is  small  compared  to 
the  peak  current  induced  in  the  horizontal  conductor  by  a  vertically 
polarized  wave  incident  at  a  small  angle  (relative  to  the  horizontal  con¬ 
ductor)  ,  The  direction  of  the  current  induced  in  the  riser  is  such  that 
a  riser  at  the  end  of  a  transmission  line  tends  to  reduce  the  maximua 
peak  current  in  the  combination  of  horizontal  and  vertical  conductors, 
but  this  reduction  would  not  be  great  for  typical- to-poor  soils. 

The  results  of  this  parametric  study  of  the  common  mode  voltages 
and  currents  induced  in  the  power  transmission  system  leave  little  doubt 
that  the  maximum  coupling  for  average- to-poor  soil  conductivity  occurs 
for  ond-on  incidence  with  vertical  polarization.  Unless  the  transmission 
lines  are  very  long  and  straight  (over  10  kilometers)  and  subjected  to  a 
vertically  polarized  wave  with  an  angle  of  incidence  of  10  degrees  or 
less,  however,  it  is  not  possible  to  achieve  the  strong  coupling  reported 
in  Baird,  Marable,  and  Nelson"^  for  conductors  of  infinite  length  %n 
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ductor  over  a  good  ground  plane.  If  the  soil  conductivity  is  10  mho/m 
or  less,  however,  the  current  induced  in  the  riser  is  small  compared  to 
the  peak  current  induced  in  the  horizontal  conductor  by  a  vertically 
polarized  wave  incident  at  a  small  angle  (relative  to  the  horizontal  con¬ 
ductor)  .  The  direction  of  the  current  induced  in  the  riser  is  such  that 
a  I'iscr  at  the  end  of  a  transmission  line  tends  to  reduce  the  maximun 
peak  current  in  the  combination  of  horizontal  and  vertical  conductors, 
but  this  reduction  would  not  be  great  for  typical- to-poor  soils. 

The  results  of  this  parametric  study  of  the  common  mode  voltages 
and  currents  induced  in  the  power  transmission  system  leave  little  doubt 
that  the  maximum  coupling  for  average- to-poor  soil  conductivity  occurs 
for  end-on  incidence  with  vertical  polarization.  Unless  the  transmission 
lines  are  very  long  and  straight  (over  10  kilometers)  and  subjected  to  a 
vertically  polarized  wave  with  an  angle  of  incidence  of  10  degrees  or 
less,  however,  it  is  not  possible  to  achieve  the  strong  coupling  reported 
in  Baird,  Marable,  and  Nelson‘s  for  conductors  of  infinite  length  in 
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practical  power  distribution  lines.  The  current  that  can  be  induced  in 

a  very  long  conductor  by  a  waye  at  grazing  incidence  has  another  linita- 

tion;  the  wave  cannot  be  maintained  near  the  surface  of  the  ground  because 

of  losses  to  the  soil  (the  powet  density  lost  to  the  soil  is  of  the  order 

T}/.]  times  that  propagating  along  the  surface) .  For  straight  transmission 
o 

lines  a  few  kilometers  long,  however,  it  is  possible  to  attain  voltages 
of  a  few  hundred  times  the  incident  field  strength  for  the  unique  end-on 
case  (^  =  p  =  0) . 

B.  Open-Circuit  Voltage  Induced  by  a  Space  yave 

The  open-circuit  voltage  induced  at  the  tcrsinals  of  a  semi-infinite 
transmission  line  is  analyzed  in  Technical  XeBprar*dum  22^  for  arbitrary 
angles  of  incidence,  soil  conductivities,  and  line  heights.  The  open- 
circuit  voltage  and  short-circuit  current  were  obtained  for  impulse  and 
step-function  incident  fields  in-  that  analysis.  Here  we  will  present  the 
transient  responses  of  the  transmission  line  for  incident  exponential 
pulses,  of  selected  decay  time  constants,  for  typical  soil  conductivities 
and  line  heights.  Although  the  high-altitude  EMP  waveform  is  usually 
represented  by  a  two-cxpor.ential  pulse,  the  exponential  rise  time  of  the 
pulse  contributes  very  little  to  the  response  of  the  transmis.  Ion  line  if 
the  rise  time-constant  is  smaller  than  the  transit  tiwe,  h/c,  between  the 
ground  and  the  tran^lssion  line.  Therefore,  the  responses  presented 
here  are  for  zero-rise- time  exponential  pulses. 

In  addition,  the  results  reported  in  Technical  Hceorandum  22'“  indi¬ 
cate  that  the  error  in  assuming  that  the  soil  behaves  as  a  good  conductor 
(i.e,,  y>  »cz)  is  very  small  for  soils  with  average  or  alx>vc  average  con¬ 
ductivity.  The  slight  error  resulting  from  this  assumption  makes  the 
open-circuit  voltages  appear  to  be  slightly  larger  than  they  would  be  if 
the  dielectric  behavior  of  the  soil  at  early  times  were  fully  accounted 
for.  Because  the  error  is  small  and  on  the  conservative  side  from  the 


E3iP  tulnorability  viewpoint,  the  results  are  given  for  the  case  when  the 
ground  is  a  good  conductor.  Neglecting  the  early-tiee  dielectric  behavior 
of  the  ground  also  siaplif ies  the  coaputation  to  the  extent  that  very 
good  results  can  be  obtained  using  tabulated  functions  and  hand  calcula¬ 
tions,  which  is  an  iaportant  advantage  when  a  quick  estiaatc  of  the  volt¬ 
age  induced  in  a  particular  transaission  line  is  required. 

The  open-circuit  voltage  developed  at  the  end  of  a  scai-infinitc 
transeission  line  of  height  h  above  a  perfectly  conducting  ground  is 

-JMt 

1  -  e  ° 

V  (o)  =  E  c  D(t,<?)  - - —  (1) 

OD  ■  O  ,  ,  t. 

jxCjtD  I-  — ) 

-t/T 

When  an  exponential  pulse  E^c  Is  incident  at  an  elevation  angle  t  and 
azi^th  angle  as  illustrated  in  Figure  1.  The  directivity  function 
D(t,9)  is  given  by 


sin  $  cos  o  (vertical 

Oc  .  .  ^  polarization) 

-i-  1  -  cos  ®  cos  o 


sin  c? 


— ^  -r  1  -  cos  ♦  cos  o 

j® 


(horizontal 

polarization) 


(2) 


and  t^  =  (2  h  sin  t)/c,  where  c  is  the  speed  of  light.  In  the  expressions 
for  D(t,9),  it  is  assuaed  that  the  phase  factor  for  the  transaissicn  line 
is  the  sase  as  that  for  free  space,  but  the  line  aay  have  sobc  attenuation, 
represented  by  cr,  if  finitely  conducting  wire  is  used  for  the  tran^ission 
line. 

The  correction  that  nust  be  applied  to  to  account  for  the  finite 
conductivity  of  the  soil  when  ?  oie  is 
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CO«»UCTOR 


FIGURE  1  COORDINATES  DEFINING  AZIMUTH  AM)  ELEVATION  ANGLES  OF 
INCIDENCE 


ilV 


E  C-Jr  DC*,o) 
o  c 


-jact 


sin  # 


( jtB  -r  — ) 


(vertical 

polarization) 


=  DC#,j?)  2  sin  ^ 


o 

c 


(horizontal 

polarization) 


Vhen  the  propagation  factor  for  the  transmission  line  is  dosinated  by 

the  capacitance  and  inductance  of  the  line  over  perfect  ground  (see 

Section  III),  the  directivity  factor  D(t,Q)  is  given  by  Eq,  (2).  In 

Eq.  (3)  above>  t  =  e  /c.  vhere  c  is  the  soil  conductivity  and  c  is 
CO'  o 

the  perwLttivity  of  free  space. 

The  total  open-circuit  voltage  is  then 


=  V 


iv 


oc 


E  c  D(# 


,v) 


-joct 


-Ject 


jo?  (j®  -T-  1/7) 


2jr  (sin  t) 

•  A 


a 


(j®  -i-  1/t) 


vherc  the  exr-onent  -i-l  is  associated  vith  horizontal  polarization  and  -1 
is  associated  *ith  vertical  polarization.  Talcing  the  Inverse  transforas, 
the  voltage  aavefora  is 


V  (t)  =  E  CTDCt,?:) 
oc  o 


(O  :£  t  5  t  ) 
o 


E  CTDC 
o 


/ 


.  \l  -t/T  .  4  (sii 

rDCf^v)  jlc  -1/  e - — 

.■/t'/T  2  ) 


(sin  f) 


c  -tVt 
e 

T 


(5) 


u 

e  dU/ 


(tS:  t  ) 

o 


*hex«  t'  =  t  -  t  .  The  first  tern  in  the  braces  is  due  to  the  geometry 
o 

(line  height  and  angle  of  incidence)  and  is  independent  of  the  properties 

of  the  soil.  TIa!  second  tern  contains  the  effect  of  the  soil  conductivity 

in  T  =6  as  veil  as  gconetric  effects, 
e  o 

The  source  iM|ccaancc  for  the  seni-iefinite  trassnission  line  is 
sinply  its  characteristic  inpcdancc  ss^Z/E.  For  typical  pover  trans- 
oission  line  confi^rations,  this  inpcdancc  deviates  very  little  front 


Z 

o 


(6) 


vhcre  a  is  the  radius  (or  effective  radius)  c£  the  pover  conductors  and 

=  /„  /c  is  the  intrinsic  inpcdancc  of  free  space.  The  Thevenin 
"0*00 

equivalent  circuit  of  the  trahssission  line  is  thus  the  open-circuit 

voltage  given  by  Eq.  (5)  in  series  vith  Z^.  The  Kortoa  equivalent  circuit 

is  a  current  source  V  /Z  in  parallel  vlth  Z  . 

oc  o  o 
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OPEN  CIRCUIT  VOLTAGE 


C,  Parametl’ic  Variation  of  the  Open-Circuit  Voltage 

The  bpen-circui?t  voltage  V  (t)  is  plotted  in  Figure  2  -for  both 

°°  -3 

polarizations  with  soil  conductivity,  varying  from  10  to  “  mho/m.  The 

waveform  ^or  o  =  “  (t  =  0)  is  produced,  by  the  first  term  of  Eq.  (5) 

e 

alone.  The  deviations  of  the  other  waveforms  from  that  for  ^  “  are 

caused  by  imperfect  reflection  of  the  incident  wave  by  the  ground.  It 

is  apparent  that  the  coupling  caused  by  the  ground  effect  is  larger  than 

the  coupling  caused- by  the  line  height  for  vertical  polarization  and  soil 

-2 

conductivities  less  than  10  mho/m.  It  is  also  apparent  that  the  ground- 
effect  coupling  to  vertical  pblarization  is  considerably  greater  than  the 


FIGURE  2  OPEN-CIRCUIT  VOLTAGE  AT  THE  END  OF  A  SEMI-INFINITE  LINE 
FOR  VARIOUS  SOIL  CONDUCTIVITIES 
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coupling  to  horizontal  polarization.  The  dir  crence  is  even  greater  than 
the  wavelorms  in  Figure  2  indicate  because  the  directivity  function 


D(\ii,cp)  is  greater  for  vertical  polarization  than  for  horizontal  polari¬ 
zation. 

The  effect  of  line  height  on  the  open-circuit  voltage  waveform  is 

illustrated  in  Figure  3^  where  the  waveforms  are  plotted  for  vertical 

-2 

polarization  incident  on  soil  with  a  conductivity  of  10  ‘  mho/m.  Also 


FIGURE  3  OPEN-CIRCUIT  VOLTAGE  AT  THE  END  OF  A  SEMI-INFINITE  LINE 
FOR  VARIOUS  LINE  HEIGHTS 


shown  as  dashed  curves  in  Figure  3  are  the  waveforms  with  <j  =  only 

this  portion  of  the  waveform  is  proportional  to  the  line  height.  For 

-2 

the  soil  conductivity,  of  -10-  mho/m^  the  ground  effect  is  larger  than 
the  height  effect,  and  changing  height  of  the  line  does  not  have  a  large 
effect  on  the  peak  open-circuit  voltage. 

The  effect  of-  incident  pulse  duration  (decay  time  constant)  on  the 
open-circuit  voltage  waveform  is  illustrated  in  Figure  4,  where  the  open- 
circuit  voltage  is  plotted  for  incident  pulse  decay  time  constants  of 
0.25,  0.5,  and  1.0  microseconds.  As  is  apparent  from  Eq.  (5)-,  the  wave¬ 
form  for  a  given  soil  conductivity,  line  height,  and  angle  of  incidence 
can  be  plotted  as  a  function  of  a  normalized  time  t/r,  in  which;  case  the 
pulse  decay  time  constant  t  affects  only  the  relative  magnitudes  of  the 
two  terms  in  braces  in  Eq.  (5).  Plotted  in  real  time  as  in  Figure  4, 
however,  it  is  apparent  that  the  pulse  duration  affects  both  the  peak 
voltage  and  the  time  required  to  reach  the  peak  value.  The  wider  the 
pulse,  the  larger  the  peak  voltage  because  the  wider  the  pulse,  the 
longer  the  segment  of  line  near  the  terminals  over  which  the  ground- 
effect  is  integrated.  Note  that  for  perfect  ground  (dashed  curves  in 
Figure  4)  the  pulse  width  has  little  effect  on  the  peak  voltage,  and  the 
open-circuit  voltage  waveform  is  essentially  the  incident  field  waveform 

except  for  the  finite  rise  time,  t  =  (2h  sin  \|f)/c. 

o 

In  Figures  2',  3,  and  4,  the  waveforms  for  vertical  polarization 

have  been  plotted  for  an  elevation  angle  ^  =  30°  and  an  azimuth  angle 

cp  =  0°,  while  the  waveforms  for  the  horizontally  polarized  incident  wave 

in  Figure  2  were  plotted  for  ^  =  30°  and  cp  =  90°  (broadside) .,  In  all 

three  illustrations,  it  was  assumed  that  or/k  «  1  and  that  of/k  was  not 

frequency-dependent.  When  Eq.  (5)  is  reexamined,  it  is  apparent  that 

changing  the  angles  of  incidence  affects  the  relative  magnitude  of  the 

±1 

ground-^effect  through  the  (sin  ^)  coefficient  as  well  as  affecting  the 


FIGURE  4  OPEN-CIRCUIT  VOLTAGE  AT  THE  END  OF  A  SEMI-INFINITE  LINE 
FOR  VARIOUS  INCIDENT  PULSE  DECAY  TIME  CONSTANTS 


directivity  iunction  D(i,cp)  and  the  delay  time  t^  for  the  arrival  of  the 

ground- reflected  wave,  in  effect,  there  are  two  directivity  functions — 

one  for  the  response  with  perfectly  conducting  ground,  D(+,(o),  and  one 

±1 

for  the  correction  term,  (sin  ♦)  The  magnitude  of  these  dir¬ 

ectivity  functions  for  vertical  polarization  are  shown  in  Figure  5  for 
various  values  of  a/k.  Note  that  for  a/k  -*  0,  both  directivity  functions 
approach  <»  as  co  0,  so  that  for  a  lossless  semi-infinite  transmission' 
line,  an  infinite  open-circuit  voltage  is  developed  for  any  finite  field 
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incident  at  cp  =  4=0.  This  situation  is,  of  course,  an  abstract  one, 
because  neither  lossless  nor  semi-infinite  transmission  lines  are  real¬ 
izable  (nor  is  the  plane  wave  to  illuminate  it).  Since  Of/k  >  0,  the 
value  of  the  directivity  functions  is  always  finite,  and  D(i,cp)/sin  t  -» 
ju)/ofc  as  ♦,  cp  0,  while  D(*,co)  I’eaches  a  maximum  value  at  some  value 
of  <1  that  is  greater  than  zero. 

The  directivity  functions  for  a  horizontally  polarized  wave  incident 
at  broadside  aro  shown  in  Figure  6.  For  ot  =  0, 


=■ 


sin  <p 


1  -  cos  ijr  cos  cp 


(7) 


which  approaches  “  as  tjf,  cp  ->  0.  For  any  real  transmission  line,  however, 
or  ^  0,  and  D(ijr,«p)  is  finite.  In  fact,  the  behavior  of  D(0,<p)  with  cp  for 


EUEVATION  ANGLE,  —  deg 

LA-79»S-10e 


FIGURE  6  VARIATION  OF  THE  DIRECTIVITY  FUNCTION  Df,(^.90J  WITH  ELEVATION 
ANGLE  FOR  or/k  «  1 


horizontal  polarization  is  identical  to  the  variation  of  for 

vertical  polarization,  since 

sin  o 

D(0,cp)  =  — - - - : -  (horizontal  polarization) 

1  -  j  -  -  cos  <p 

K 

(8) 

sin  6 

DC^jO)  = - — - ^ -  (vertical  polarization) 

1  -  j  ^  -  cos 

Thus  the  dashed  curves  of  Figure  5  may  be  used  to  illustrate  the  behavior 
of  D(0,<p)  for  horizontal  polarization  if  f  in  Figure  5  is  replaced  by  cp. 
Figure  6  illustrates  that  the  directivity  function  for  the  ground  effect, 
sin  if  DCv,0),  illustrated  by  the  solid  curve,  is  always  equal  to  or  less 
than  the  dii'cctivity  function  D(^t,0)  for  the  wire  over  perfect  ground. 

Uccausc  the  open-circuit  voltage  is  dominated  by  the  ground  effect 
for  vertical  polarization  and  average  or  Icss-than-averagc  soil  conduc¬ 
tivity,  and  because  the  response  to  vertical  polarization  is  greater  than 
the  response  to  horizontal  polarization,  the  maximum  open-circuit  voltage 
will  be  developed  when  the  transmission  line  is  subjected  to  a  vertically 
polarized,  end-on  illumination  (i?  =  cp  =  0),  Then  D(t,<p)/sin  jk/(Y  -  jk); 
and  the  problem  becomes  one  of  a  vertically  polarized  wave  propagating 
along  the  transmission  line.  It  is  no  longer  reasonable,  iti  this  situa¬ 
tion,  to  assume  a  unifora  plane  wave  and  a  semi-infinite  line  because 
thei-e  exists  no  practical  case  in  which  a  uniform  field  can  be  generated 
over  a  very  long  line  when  the  wave  is  incident  at  a  grazing  angle  (f  =  0) . 
For  the  case  where  f  =  9  =  0,  therefore,  we  will  assume  that  the  total 
vertical  electric  field  varies  as  1/z  and  that  the  length  of  the  trans¬ 
mission  line  is  finite. 


D.  Vertically  Polarized  Wave  at  Grazing  Incidence 
1.  General 

A  case  for  which  a  vertically  polarized  wave  incident  at  ^  =  cp 
=  0  is  of  interest  is  illustrated  in  Figure  7.  In  this  case^  the  BMP  is 
generated  by  a  surface  burst  on  the  transmission  line  (or  an  extension 
of  the  line) .  Near  the  source  point  (in  the  close-in  region) ,  the  radia¬ 
tion  from  the  weapon  will  ionize  the  air  sufficiently  to  cause  it  to 

behave  as  a  conductor  (a^  uie  ),  and  further  away  (in  the  intermediate 

o 

region),  the  ionization  will  be  weaker  but  still  sufficient  to  produce 
significant  attenuation  of  signals  propagating  away  from  the  source  on 
the  line.  Although  it  is  evident  that  there  will  be  coupliiig  to  the 
transmission  line  in  the  close-in  and  intermediate  regions,  the  analysis 
of  the  coupling  in  these  regions  is  beyond  the  scope  of  this  report. 

Here  we  discuss  only  the  coupling  that  occurs  outside  these  regions, 
where  the  BMP  wave  behaves  as  a  vertically  polarized  wave  propagating 
outward  from  the  source  at  the  speed  of  light.  The  analysis  is  performed 


FIGURE  7  PHYSICAL  ARRANGEMENT  THAT  PRODUCES  A  VERTICALLY  POLARIZED 
WAVE  WITH  END-ON  INCIDENCE 
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in  such  a  manner  that  the  signals  induced  in  that  part  of  the  line  closer 
to  the  source  can  be  superimposed  on  the  results  calculated  for  the  region 


Jurther  from  the  source. 

For  our  analysis^  the  transmission  line  axis  is  assumed  to  pass 
through  the  source  point  at  z  =  0  and  is  assumed  to  leave  the  intermediate 
region  at  z  and  extend  to  the  terminals  at  z  >  z  as  illustrated  in 

x  A  X 

Figure  7.  The  wave  emanating  from  the  intermediate  region  is  assumed  to 

propagate  as  though  the  ground  were  a  perfect  conductor  so  that  E  =  T|  H  , 

*  p-  y 

where  T)  is  the  intrinsic  impedance  of  free  space,  and  E  =  E  e~j^^/z, 
o  •  '  X  o 

To  determine  the  z-component  of  the  field  in  the  ground,  however,  we  will 

concede  that  the  ground  is  not  a  perfect  conductor  and  observe  that  for 

a  finitely  conducting  ground,  the  z-componeht  of  the  electric  Tield  at 

the  surface  will  be  K  (0,z)  =1)11  (0,z),  where  Tl  is  the  intrinsic  imped- 

z  y 

anco  of  the  soil,  and  T|  <3C  Since  we  will  need  the  z-component  of 

the  electric  field  at  the  wire  height  to  calculate  the  current  induced 
in  the  wire,  we  can  approximate  this  field  by 

(9) 

X  =  0 

where  5E  /f>x  is  obtained  from  the  curl  E  at  the  surface: 
z 


5e 

E  (h,z)  »E  (0,z)  +  h 
z  ^  z  ox 


dE 

2 

5x 


=  jtuii  H 
o  y 


5E 
^  X 

dz 


JUJU 

o 


E  + 

X 


(-jk  -  -) 
z 


E 

X 


(10) 
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since  E  =  E  and  H  =  E  /Tl  .  The  field  at  the  wire  height  is 

X  o  y  X  o 

thus 


(1^) 


In  dei’iving  this  field,  we  have  (1)  neglected  attenuati^^n  of  the  wave 
due  io  losses  to  the  ground,  and  (2)  applied  a  first  order  correction 
to  the  field  at  the  surrface  to  obtain  the  field  at  the  wire  height.  In 
addition  it  should  be  pointed  out  that  E^  is  now  the  total  vertical  elec¬ 
tric  field,  so  that  it  corresponds  to  2E^  in  the  space  wave  case  discussed 
in  Section  II-B. 

2.  Current  Induced  in  Wire 

The  current  induced  in  the  wire  aay  now  be  obtained  Tor  the 

region  between  the  edge  of  the  internediate  region  at  and  the  tcminals 

at  z  .  This  current  is  obtained  Iron  the  differential  equation 
2 

2 

^  =  -  y  E^(h,z)  (12) 

dz*”  * 


which  has  a  solulion 

I(z)  =  [K  •I'  P(z)3  e~^^  -i-  Ck  -i-  Q(z)]  e^^ 

X.  M 


where 


o 


E  (h,v)  dv 
z  ' 


(13) 
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Q(z) 


o  z 


E  (h.v)  dv 
z 


K.  = 


-yz  yz 

2  2 
-  Q(Zj^)  e 


v(vV 


—  0  D  C 


-Y(z  -z  ) 
2  i 


K_  = 


^>2  ® 


YZ^  -YZ^ 

-YZg  Pj^QCZj^)  e  -  PCZg)  e 


Y<Zj,-Zj) 


•"  D  D  © 


-Y(z2-z^) 


Z.  -  Z 
X  o 


"  Z.  Hr  Z 


y  =  Of  +  jp 


We  will  assume  that  the  coupling  to  the  line  extending  into  the  inter¬ 


mediate  region  can  be  represented  by  a  current  source  and  an  impedance 


Z  at  z  =  z  as  illustrated  in  Figure  8,  The  effect  of  this  source  can 
o  1 


be  superimposed  on  the  induced  current  calculated  here  for  that  portion 


of  the  line  between  z^  and  outside  the  intermediate  region.  The  short 


circuit  current  at  the  terminals  at  z^  due  to  the  coupling  between  z^ 


and  z  will  be  obtained  by  letting  the  impedance  at  z  be  zero,  so  that 
2  2 


p  =  -1.  Because  the  source  impedance  at  z,  is  Z  ,  p  =0  and  K  =  0, 
2  1  o'  *^1  1 


Also,  because  the  short  circuit  current  at  z^  is 


Kz^)  =  o 


-Yz^  Y=^2 

+  K^c 


(14) 
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FIGURE  8  GEOMETRY  FOR  ANALYSIS  OF  THE  COUPLING  CONFIGURATION  OF 
FIGURE  7 


where 


K  =  e  P(z  ) 

2  2 


or 


-Y2„ 


Kz^)  =  ^PCz^)  e 


When  E  (h.z)  is  substituted  into  the  integral  for  P(z),  the  short  circuit 
z 

at  the  terminals  induced  between  z  and  by  the  wave  is 

1  A** 


*  o  ^ 


E  e 

V  =  ^ 


^  f 

T)  J 

o  z 


2  (Y-jk)v 
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dv  -  h 


(Y-jJ«)v 


dv 


E  h  e'^^2  (  ^Y-Jk)z^ 

o  }e  e 


—  +  (iif  - 

1  L  o  J 
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2  (Y'*jk)v 
X  /  - - -  dv! 


(15) 
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The  coupling  to  the  wire  can  now  be  resolved  into  three  com¬ 
ponents  according  to  the  nature  of  coupling.  If  the  wire  and  soil  were 
perfect  conductors,  then  I)  =  0  and  Y  =  and  the  short-circuit  current 


in  the  wire  at  z  would  be 
2 


I 

B 


(16) 


The  coupling  in  this  case  is  caused  by  the  nonuniformity  in  the  magni¬ 
tude  of  the  incident  electric  field  along  the  wire,  which  causes  the 
potential  at  the  wire  height  at  z^  to  differ  from  the  potential  at  the 
wire  height  at  z^. 

If  Y  =  and  T]  r  0,  we  obtain  an  additional  component. 


E  T|  -jk*  ^2 

00  1 

v/hich  is  caused  by  the  finite  conductivity  of  the  soil  (1)  yju>p/o) . 

This  component  is  the  integral  (between  z  and  z  )  of  the  z-component  of 

X  2 

the  electric  field  induced  in  the  finitely  conducting  soil  by  the  inci¬ 
dent  wave.  If  exists  in  some  form  whether  or  not  the  incident  wave  has 

geometric  (1/z)  variation,  although  the  log  (z  /z  )  term  is  character- 

2  1 

istic  of  the  1/z  variation  of  the  incident  field. 

Finally,  if  Y  /  third  component  is  obtained.  This  com¬ 

ponent,  caused  by  the  disparity  in  the  propagation  characteristics  of  the 
incident  wave  and  the  induced  responses,  modifies  both  the  geometric 
component  and  the  gi’ound-fiold  component  (since  y  r  jk)  and  adds  to  these 
a  component. 


I 


^CY-jk)v 

w 


dy 


(18) 


which  exists  only  if  there  is  a  disparity  betweea  the  phase  factors  of 
the  incident  wave  and  those  of  the  tran^ission  line.  This  phase-dispari'fcy 
coeppnent  is.  in  reality^  a  gcoeetric  coeponcat  also^  since  the  integral 
in  Eq.  (15)  whose  coefficient  is  y-jk  occurs  because  of  the  l/z  varia¬ 
tion  of  the;  incident  fxeld.  ' 

For  typical  soil  ccnductivitfos  and  line  lengths,  the  gcoaetric 
term  given  by  Eq.  (16)  is  significant  only  at  very  low  frequencies  unless 
the  line  is  very  high  or  very  close  to  the  source,  and  the  phase-discrepancy 
tera  given  by-  Eq.  (18)  is  ssall  <x>aparcd  to  the  ground-effect  texs  given 
by  Eq.  (17).  The  short-circuit  current  induced  in  the  transaission  line 
by  an  exponential  pulse  E  /z  c  '  is  tfuis 


when  phase  is  referred  to  the  tereinals  at  z^. 
The  current  wavefora  is  thus 
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NORMALIZED  CURRENT 


for  a  zero-rise-time  exponential  pulse.  For  a  two-expohential  pulse  with 
non-zero  rise  time,  the  current  waveform  can  be  obtained  by  superposing 
two  solutions  such  as  in  Eq.  (20) .  The  waveforms  for  single  and  double 
exponential  transient  fields  are  shown  In  Figure  9  for  a  pulse  decay  time 
constant  of  1  ps  and  soil  conductivities  varying  from  10"  to  10"  mho/m. 
For  the  two-exponential  pulse  (solid- .cui’ves)_,  the  rise  time  constant,  is 
4  ns.  The  early  time  and  late  time  approximations  do  not  meet  at  t  = 
for  the  lower  conductivities,  so  they  have  been  joined  by  a  vertical 
line.  The  dashed  curves  are  the  zero  rise  time  (single  exponential) 
responses. 


NORMALIZED  TIME,  Mr 


LA-7995-1 1 1 


FIGURE  9  SHORT-CIRCUIT  CURRENT  AT  THE  TERMINALS  OF  A  LINE  1000  METERS 
LONG  AND  1000  METERS  FROM;  THE  SOURCE  WITH  1  V/m  INCIDENT 
ON  THE  END  NEAREST  THE  SOURCE  (Z^  =  300  =  1000  m, 

Z2  =  2000  lii) 


31 


The  waveforms  of  Figure  9  are  such  short,  high-amplitude  pulses 
that  it  was  necessary  to  plot  them  on  a  log-log  scale  to  view  significant 
details.  It  is  noted  that  the  peak  value  and  the  tine  to  reach  the  peak 
value  both  increase  as  soil  conductivity  decreases.  This  is  because  the 
poorer  conducting  soil  behaves  as  a  dielectric  longer  and  supports  a 
larger  z-component  of  electric  field  when  it  behaves  as  a  conductor. 

Also  interesting  is  the  fact  that  the  short-circuit  current  waveform  has 
a  negative  undershoot  beginning  at  about  one  decay  time- constant. 

For  comparison  of  the  end-on  illumination  of  the  finite-length 

transmission  line  with  the  oblique  incidence  on  the  semi-infinite  line, 

we  will  compute  the  peak  current  and  time-to-peak  for  a  soil  conductivity 
-2 

of  10  mho/m  and  an  exponential  pulse  decay  time  of  1  ps.  Assume  the 

Tinite-length  line  is  1  kilometer  from  the  source  and  1  kilometer  long, 

and  that  the  field  strength  at  the  end  nearest  the  source  is  -1  V/m 

(E  =  1000),  The  characteristic  impedance  of  the  line  is.  taken  to  be 
O; 

300  ci-ns.  Then  the  peak  short  circuit  current  from  Eq,  (20)  and  Figure 
9  is  0,65  A  occurring  at  t  10”®  seconds.-  From  Figures  4  and  5,  the  loss¬ 
less  semi-infinite  transmission  line  illuminated  with  a  plane  wave  of 

1  V/m  at  f  =  10°  will  have  a  peak  short-circuit  current  of  about  2.C  A 

-7 

occurring  at  t  »  6  X  10  seconds. 

Although  some  of  the  difference  between  these  two  results  can 
be  attributed  to  the  non-uniform  illumination,  most  of  it  can  be  attrib¬ 
uted  to  the  finite  length  of  the  line  in  one  case  and  its  infinite  length 

in  the  other.  On  the  semi-infinite  line,  the  peak  current  occurring  at 
-7 

t  =  6  X  10  seconds  contains  current  that  was  induced  ct/(l  -  cos  ^) 

=11,8  kilometers  away  from  the  terminals.  Thus,  unless  the  traiismission 
line  is.. straight  and  at  least  11.8  kilometers  long,  the  short-circuit 
current  will  never  attain  the  peak  value  predicted  for  a  semi-infinite 
line.  One- should,  therefore,  use  caution  in  interpreting  ishe  currents 
and  voltages  predicted  for  semi-infinite  lines  with  small  angles  of  inci- 
ence,'^ 
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E.  Response  of  Vertical  Element 

Vertical  el’ements  such  as  ground;  wires  and  service  entrance  conduits 
also  interact  with  the  incident  wave  and,  thus,  have  current  induced  in 
them.  Because  only  the  vertically  polarized  wave  has  a  component  of 
electric  field  in  the  vertical  direction,  the  vertical  elements  actively 
interact  with  only  the  vertically  polarized  wave;  for  horizontally  polar¬ 
ized  waves,  the  vertical  .elements  behave  as  passive  impedances  with  delay 
times  associated  with  their  length.  The  vertical  element  is  considered 
to  be  a  biconic  transmission  line  with  its  upper  end  terminated  in  its 
characteristic  impedance  and  its  lower  end  short-circuited  to  the  ground. 
The  current  induced  in  this  biconic  transmission  line  by  the  vertical 
component  of  the  electric  field  can  then  be  determined  in  the  same  way 
as  the  current  in  the  horizontal  transmission  line  over  ground.  As  was 
pointed  out  in  Technical  Memorandum  22,^  the  current  induced  in  the  verti¬ 
cal  riser  opposes  the  current  induced  in  the  horizontal  conductor  when 
|fp|  <  Tr/2  and  aids  it  when  |<p|'.  >  it/2.  Since  the  coupling  to  the  hori¬ 
zontal  conductor  is  greatest  for  |q)|  <  Tr/2 — when  D(i|(,(p)  is  large — the 
current  induced  in  the  vertical  element  tends  to  reduce  the  -maximum  cur¬ 
rent  calculated  for  the  horizontal  conduction  alone. 


For  a 
a  vertical 


vertically  polarized  wave  incident  at  an  elevation  angle  ijr  on 
element  of  height  h,  the  current  at  the  top  will  be 


1(h)  = 


E  c  cos  ill  , 
o  1 


1  -  e 


.  2h  . 

-Jiu  —  (1  +  sin  ijf) 
c 


2Z  jw  L  1  +  sin  \jr 


(21) 


2h  .  ^  2h 

-ju)  sin  y  -j(ju  — 
c  c 

0  -  0 


1  -  sin  f  J  juj  -I-  1/t 


I 


3o 


for  an  exponential  pulse  E  e  .  In;  the  expression- above.  Z  is  the 

o  o 

characteristic  impedance  of  the  ibiconic  'transmission  line  given  by 


Z 

o 


•o  ,  2h 
»  —  log  — 

2n  a 


(22) 


and  the  ground  is  assumed  to  be  a  perfect  conductor  (R  =  1).,  The  wave- 

V 

form  of  the  current  at  the  top  of  the  vertical  element  is 


i(h,t) 


whore  t  ^  2h/c  (1  +  sin  i)).,  t  =  2h/c  sin  and  t  =  2h/c  . 

1  2-  3'  ^ 

For  a  vertically  polartzed^wave  at  grazing  incidence  (i)  =  0), 


i(h,t) 


E  CT 
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rt/T 
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-(t^t  )/t\ 
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3 


(24) 


Plots  of  the  current  waveform  are  shovm  in  Figure  10  for  ^  =  0  and 
=  30®  for  an  incident  exponential  pulse  with  a  decay  time  constant 
T  =  1  ps  and  a  riser  height  h  =  10  m.  The  waveforms  for  these  two  angles 
of  incidence  are  quite  similar,  but  the  leading  edge  of  the  waveform  for 
i;  =  30®  is  more  complex  because  of  the  difference  in  time  of  arrival  of 
the  upward  and  downward  traveling  current  waves  induced  by  the  direct  and 
ground-reflected  electric- field  waves. 


A  comparison  of  the  current  induced  in  the  riser  for  conditions 
similar  to  those  used  earlier  for  comparing  the  current  between  the 
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FIGURE  10  CURRENT  DELIVERED  TO  A  MATCHED  LOAD  AT  THE  TOP  OF  A 
VERTICAL  rIsER  BY  A  VERTICALLY  POLARiZED  INCIDENT  WAVE 

horizontal  -conductor  of  finite  length  and  the  Bmi-infinite  conductor 

shows  that  the  peak  current  induced  in  the  riser  by  a  1- V/m  Incident 

exponential  pulse  (i|r  =s  0,  h  =  10  m,  T  =  1  ps)  is  about  0.064  A.  The  peak 

current  in  the  finite-length  horizontal  conductor  was  0.65  A  tuid,  in  the 

-2 

semi-infinite  conductor,  it  was  2,8  A  (a  =  10  mho/m) .  The  peak  current 
induced  in  the  riser  is  thus  only  about  10  percent  of  the  peak  ground- 
effect  current  induced  in  a  1000-meter-lohg  horizontal  conductor,  but  it 
is  flowing  in  the  opposite  direction  from  the  ground-induced  current. 

For  perfectly  Conducting  ground,  however,,  the  peak  current  induced  in 
the  riser  is  one-half  as  large  as  that  induced  in  a  semi-infinite  hori¬ 
zontal  conductor,  so  that  the  higher  the  soil  conductivity,  the  more 
significant  the  riser  current. 


SECTION  III 

PROPAGATION  OF  CURRENTS  INDUCED  IN  POWER  LINTIS 

A.  Introduction 

In  previous  reports  tha  coupling  of  incident  electromagnetic  energy 
onto  power  lines  has  been  analyzed  and  the  analysis  confirmed  with  mea¬ 
surements.  This  report  discusses  how  these  induced  currents  propa¬ 

gate  down  a  power  line  over  real  earth  and  with  discontinuities.  The 
induced  current  will  be  changed  by  attenuation  and  dispersion  due  to  the 
real  earth.  In  addition,  the  waveform  will  be  altered  by  a  variety  of 
discontinuities  that  occur  in  actual  power  lines — bends,  junctions,  spurs, 
and  grounded  neutrals.  Since  propagation  is  important  for  wavelengths 
that  are  short  compared  to  the  height  above  ground,  the.  transmission  line 
solutions  are  not  strictly  correct.  In  this  chapter,  the  significance 
of  these  effects  is  experimentally  evaluated.  In  the  final  section  some 
consideration  is  given  to  coupling  -to  multi-wire  lines. 

B,  Real  Earth  Effects 

The  propagation  of  a  signal  down  a  wire  over  real  earth  has  been 
worked  out  by  Sunde.^  An  analysis  that  gives  similar  results  was  pre¬ 
sented  in  Technical  Report  5.^  We  shall  not  repeat  these  analyses  but 
will  give  the  results  from  Sunde's  analysis. 

The  propagation  constant,  y,  is  given  by 

Y  =  Y  H(ju)) 
o 

where  Y  =  jk  =  j(u)/c)  free-space  propagation  constant,  c  =  velocity  of 
o 

light  in  vacuum,  and  H(ju))  is  a  function  that  contains  the  earth  and  the 
line  parameters: 
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H(j(u)  = 


1  +  {&n  Mn 


1  + 


1/2 


/•  ,  .1/2,. 


1/2 


(25) 


For  most  cases  of  interest  the  second  term  on  the  right  is  small 
compared  to  1  so  that 


H( ju))  s:  1  +  1/2 


1  + 


1/2 


.1/2 


1/2 


(26) 


T^=uah 


T  =  pa  a 
a  i 


where  h  =  line  height  over  ground,  a  =  line  radius,  p  =  permeability. 


a  =  soil  conductivity,  and  a  ~  wire  conductivity.  The  terms  JjuiT  and 

i  h 

yjtUT  may  be  recognized  as  complex  and  related  to  the  line  height,  h,  and 
conductor  radius,  a,  normalized  to  the  appropriate  skin  depth,  since 

yjujT^  =  -^(l  +  jl)  6  =  soil  skin  depth  = 


u)pa.\-l/2 


yjujT  =  —  (1  +  jl)  6  =  wire  skin  depth  =  1 - 1 

a  0.  i  \  2  / 


To  gain  some  insight  into  the  propagation  over  real  earth,  Sunde 
calculates  the  step  function  response  of  a  line.  Although  there  is  no 
closed  form  solution  for  all  times,  an  analytic  solution  may  be  found 
for  both  small  times  and  for  large  times,  where  the  distinction  between 
small  and  large  times  depends  upon  the  value  of  h/6.  If  h/6  »  1,  the 
small-time  solution  is  appropriate.  If  h/6  «  1,  the  large-time  solution 
is  appropriate. 
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For  small  times,  the  current  at  a  point  z  down  the  line  normalized 


to  the  current  at  z  =  0  is 


(z,jaj)  jojz 

— =  exp - — 

,0<u)  c 


where 


(27) 


For  most  cases,  t  »  ^  so  that  t  ««  t  .  The  transform  of  this 
ah  1  h 

expression  for  current  when  a  step  function  of  current  is  introduced  at 
z  =  0  is 


I(2,t) 

I(o,tJ 


erfc 


(28) 


for  small  times,  where  erfc  (z)  =  complimentary  error  function  of  z,  and 

t  =  t  -  z/c. 
o 

For  large  times  the  normalized  current  in  the  frequency  domain  is 


^  exD  -  ^ 
I(o,ju))  c 


JUJZ 


.  ^  2h 
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The  time  domain  solution  with  step- function  input  is 


I(z,t) 

I(o,tJ 


Vl/2 


1  +  2 


vnr 


-  2h 
4c  t  in  — 
o  a 


(30) 


for  large  times. 
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For  many  cases  2(t  /nr  )  1.  which  further  simplifies  the 

o  c. 

solution.  Note  that  for  large  times  the  response  is  independent  of  the 
value  of  soil  conductivity  and  only  weakly  dependent  on  the  height  above 
ground  through  the  logarithmic  term  in  h/a.  This  contrasts  with  the 
small  time  solution,  which  is  a  function  of  both  o  and  h. 

The  contrast  between  small-  and  large-time  solutions  can  be  illus¬ 
trated  by  looking  at  the  attenuation  constant  as  a  function  of  frequency 
for  different  values  of  cf  and  h  as  shown  in  Figures  11  and  12.  This  con¬ 
stant  is  derived  from  the  expression  for  K( jco) .  The  values  for  which 
computations  were  carried  out  encompassed  the  parameters  for  which  mea¬ 
surements,  which  are  described  in  a  subsequent  section,  were  made.  These 
same  charts  can  be  used  for  other  ranges  of  parameters  by  employing  the 
following  scaling  relations.  For  a  line  height  and  wire  radius  k  times 
greater  than  those  given  in  Figures  11  and  12,  divide  the  frequency,  atten¬ 
uation  and  conductivity  scales  by  k.  If  the  wire  radius  is  not  scaled 
exactly,  the  attenuation  constant  should  also  be  multiplied  by 

(In  2h/a)  .  /In  2kh/a^  where  a.'  is  the  actual  radius.  The  values  for 
figure 

100:1  scaling  (full  scale  values)  are  shown  in  the  parentheses. 

In  Figures  11  and  12,  the  attenuation  due  to  the  wire  has  been 
omitted  since  it  changes  the  result  by  only  a  few  percent. 

The  large-time  solution  is  most  prominent  in  Figure  11  at  frequencies 
below  about  100  The  attenuation  constant  varies  less  than  a  factor 

of  2  when  is  increased  100  times.  The  variation  of  ct  with  frequency 
is  almost  linear. 


In  Figure  12,  where  the  height  is  increased  by  a  factor  of  6,  the 

attenuation  constant  depends  much  more  strongly  on  't,  approaching  a 
-1/2 

<7  dependence  at  very  high  frequencies.  In  this  case,  or  also  varies 
1/2  1 

with  f  ,  not  with  f  as  previously  mentioned. 


39 


FREQUENCY  —  Hz 


(10^)  no®)  no®)  no’)  <io®) 


FIGURE  n  ATTENUATION  CONSTANT  OF  A  LINE  OVER  GROUND  AS  A  FUNCTION  OF 
FREQUENCY,  WHERE  h  =  17  FEET  FULL  SCALE 

A  conparison  of  Figures  11  and  12  also  shows  that  at  low  frequencies 
a  is  relatively  independent  of  h,  while  at  high  frequencies  a  varies 
approximately  as  h 

If  the  complete  expression  for  y  is  manipulated  to  extract  the  atten 
uation  and  phase  constant  separately,  rather  simple  expressions  result 
when  the  wire  loss  is  neglected. 


tu/iMdiu 


(U/tMdau 


These  expressions  are  accurate  for  the  entire  frequency  range  of 
applicability  of  Equation  26.  The  key  parameter  is  seen  to  be  h/b,  since 
In  2h/a  is  a  slowly  varying  quantity  with  line  parameters.  The  only  other 
parameter  is  k  =  U)/c,  which  is  sixply  rela.ted  to  the  frequency. 

The  asymntotic  behavior  of  a  when  h/6  is  very  small  and  very  large 
compared  to  1  illustrates  how  the  line  and  soil  parameters  affect  the 
small-  and  large-time  solutions.  When  h/6  1,  the  arc  tangent  term 

reachet  a  constant  value  of  n/4.  Therefore 


O' 


kn 


^  2h 
8  in  — 
a 


(33) 


for  late  times.  For  this  condition  o  is  independent  of  a  and  only  weakly 
dependent  on  h  through  the  In  2h/a  term;  or  is  linearly  dependent  on  fre¬ 
quency. 

When  h/6  »  1,  the  arc  tangent  term  reaches  l/(2h/6)  so  that 


O' 


4h/6  in  ~ 
a 


(34) 


for  early  times.  Under  these  conditions  a  varies  with  a  and  with 
1/2 

f  .  In  addition,  a  is  inversely  proportional  to  h  when  the  variation 
of  in  2h/a  is  neglected. 

A  good  analytic  approximation  for  all  frequencies  is 


8.2  X  lo"^  f 


in 


2h 


[1  +  1.47  h/6] 


a  = 


nepers/m 


(35) 


C-.  Measurements  of  Pi'opagation  Effects 

To  confirm  the  above  calculations,  especially  for  wavelengths  com¬ 
parable  to  the  line  height,  two  types  of  measurement  were  made.  The  step 
function  response  of  a  line  over  real  earth  was  measured  and  compared  to 
the  calculated  values.  This  comparison  was  effective  in  checking  per- 
fonnances  at  the  higher  frequencies.  Because,  however,  the  attenuation 
falls  to  such  low  values  at  low  frequencies,  this  method  does  not  insure 
detection  of  an  error  in  the  low  frequency  value  of  a.  In  addition, 
therefore,  Q  was  measured  with  CW  signals  over  the  lower  range  of  fre¬ 
quencies. 

D.  Step  Function  Response 

In  the  experimental  set-up,  which  is  shown  in  Figure  13,  a  mercury 
reed  pulse  generator  produced  an  approximate  step  function  with  a  rise 
time  of  about  300  ps.  This  was  fed  into  a  balun,  which  in  turn  excited 
a  transmission  line  above  ground.  The  line  was  generally  terminated  in 

but  since  the  significant  measuring  time  was  shorter  than  the  time 
required  for  end  reflections  to  arrive  at  the  measuring  point,  the  ter¬ 
mination  point  did  not  affect  the  results.  A  CT-1  current  transformer 
was  used  to  measure  the  current  waveform  along  the  wire. 

The  earth  was  salted  to  increase  its  conductivity  to  between  0.2  and 
0.3  mho/m  as  measured  by  a  test  cell  at  1  MHz. 

The  measured  and  calculated  results  for  a  single  wire  line  2  Inches 
above  ground  and  16  feet  from  the  feed  point  are  shown  in  Figure  Id. 

The  calculated  value  is  for  the  late  time  solution.  The  agreement  Is 
quite  good.  The  early  time  solution  gives  quite  a  different  waveform, 
which  differs  considerably  from  the  measured  waveform.  For  our  waveform 
most  of  the  energy  is  below  1  GHz.  At  this  frequency  h/6  1.  Thus,  for 

most  of  the  frequencies  h/6  <  1,  and  the  late  time  solution  should  apply. 
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FIGURE  13  EQUIPMEMTAL  SET-UP  FOR  MEASUREMENT  OF  CURRENT  Y/AVEFORMS 
OF  A  LINE  OVER  GROUND 


In  jontrast,  the  sa=e  scasureccnts  Esade  on  a  single  s-ire  line  at  a 
he.'.ght  of  12  inches  sho«ed  hardly  any  degradation  in  waveshape  after 
t  avellag  16  feet  froa  the  feed  point.  For  this  configuration  h/fi  >  1 
fcr  fr^^quencies  down  to  10^  Hz,  and  the  snail  tice  solution  was  appro¬ 
priate.  This  solution  showed  no  significant  degradation  in  waveshape  at 
this  distance.  At  a  distance  of  40  feet,  snail  changes  could  be  neasured 
that  were  in  accord  with  the  snail  tine  solution  (sec  Figure  15). 

Sioilar  ceasurenents  were  nade  with  three-wire  lines  at  both  2-inch 
and  12-inch  heights.  The  results  are  shown  in  Figures  16  and  17.  The 
agreenent  between  the  ccasured  and  calculated  wavefoms  is  quite  good 
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FIGURE  14  CURRENT  WAVEFORM  OF  A  SINGLE  WIRE  LINE  1600  FEET  (full  scale) 
FROM  THE  FEED  POINT,  WHERE  h  =  17  FEET  FULL  SCALE 


when  the  small  Lime  solution  is  used  for  h/6  >  1  and  the  lurKC  time  solu¬ 
tion  when  h/6  <  1. 

The  three-wire  solution  is  essentially  the  same  as  the  single  wire 
solution  with  the  exception  that  the  wire  radius  I'or  the  multi-v/iro  line 
is  an  effective  radius  that  is  larger  than  the  individual  v/ire  radius 
and  smaller  than  the  total  distance  between  the  wires.  This  difiorence 
in  the  effective  radii  accounts  for  the  smaller  impedance  of  the  three- 
wire  lino  over  ground. 

The  conclusion  from  these  measurements  is  that  Sunde's  formulation 

is  accurate  for  predicting  the  change  in  waveshape  of  a  signal  down  a 

transmission  line.  Since  the  formulation  does  not  include  any  radiation 

loss  terms,  it  is  concluded  that  any  radiation  loss  is  negligibly  small 

compared  to  t.  loss  due  to  the  finite  conductivity  of  the  ground.  Our 

measurements  at  a  height  of  12  inches  corresponded  to  a  height  of  one 

wavelength  at  a  frequency  of  1  GHz.  These  measurements  may  be  considered 

a  20:1  scale  model  of  a  three-wire  line  20  feet  high  over  soil  with  con- 

_2 

ductivity  of  about  10  mho/m;  the  conclusion  that  radiation  losses  for 
pulses  with  rise  times  of  greater  than  6  ns  are  negligible  still  holds. 
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Because  the  attenuation  at  low  frequencies  is  too  low  to  be  detect¬ 
able  in  step  function  measurements,  the  attenuation  at  these  frequencies 


O'  w  +  p 


"i” 


(36) 


where  w  =  distance  from  short  circuit,  =  absorption  term  o. 

for  a  perfect  short),  P  =  2n/\,  and  Lv  =  3  dB  width  of  currem 

This  technique  allows  very  small  values  of  attenuation  to  be  i 

Large  plates  were  used  as  a  short  circuit.  These  kept  the  va‘. 

low  by  providing  a  low  impedance  to  ground  and  by  minimizing  i 

from  the  end  of  the  line.  Measurements  were  made  over  a  freqt 

of  10  to  100  MHz  for  a  line  2  inches  above  ground,  which,  on 

ratio,  corresponds  to  a  full-scale  line  at  a  height  of  20  fee 

the  full-scale  frequencies  correspond  to  83  to  830  kHz.  The  i 

-3 

conductivity  corresponds  to  about  2  X  10  mho/m. 

The  measured  results  are  shown  in  Figure  11,  which  also  j 
culations  based  on  Sunde's  analysis.  Wire  conductor  losses  ai 
losses  are  assumed  negligible  in  the  calculation,  A  perfect  f 
cuit  was  assumed  in  reducing  the  measurements.  At  the  higher 
it  was  possible  to  measure  the  null  width  at  more  than  one  nui 
way  was  estimated  and  found  to  be  small  at  these  frequenci' 

When  a  conductivity  of  about  0.2  mho/m  is  assumed,  it  cai 
that  the  measured  values  are  slightly  higher  than  the  calcula 
(at  most  30  percent  higher)  and  follov/  the  theoretical  variaL 
frequency  quite  well. 

The  phase  constant  p  was  also  checked  by  measuring  the  pi 

the  null.  With  a  perfect  short  circuit  and  free  space  propag; 

current  nulls  should  occur  at  (2nfl))v  /4.  The  deviation  from 

o 

tions  is  proportional  to  the  deviation  in  phase  constant.  Fi{. 
shows  the  measured  and  calculated  values  of  the  deviation  froi 
phase  constant  for  the  line  2  inches  above  ground.  The  point.' 
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reasonably  close  to  the  value  for  a  conductivity  0.2  mlio/m.  The  measured 
value  of  conductivity  was  between  0.2  and  0.3  niho/m. 


SECTION  IV 


DISCONTINUITIES  ON  POWER  LINES 

The  signals  coupled  to  the  power  lines  will  not  always  pi 
straight  to  the  load.  There  may  be  a  variety  of  discontinuiti 
ing  bends,  junctions,  and  spurs.  In  this  section,  we  present 
retical  and  experimental  results  for  propagation  of  signals  wl 
tinuities  are  present. 

Analysis  is  usually  restricted  to  the  case  where  the  lim 
small  compared  to  the  wavelength.  Thus,  one  purpose  of  the  mt 
was  to  see  if  significant  deviations  from  theory  could  be  four 
line  height  was  comparable  to  the  v/avelength. 

A.  Bends 

The  simplest  discontinuity  is  a  bend  in  the  line.  King, 
this  type  of  discontinuity  when  gh  <3C  1,  considered  the  geomei 
in  Figure  20 J  His  analysis  shows  that  the  bend  has  the  effet 
ing  the  inductance  Ji(w)  and  capacitance  c(w)  per  unit  length  < 
the  region  around  the  bend.  The  values  of  -6(w)  and  c(w)  are 
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FIGURE  20  TWO-WIRE  LINE  WITH  A  BEND 
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Measurements  were  made  with  a  single  wire  line  1  foot  and  2  feet 
above  ground  for  a  range  of  bend  angles  from  0°  to  157,5°.  A  step  func¬ 
tion  signal  with  rise  time  of  about  300  ps  was  fed  onto  the  line.  A 
current  transformer  was  mounted  1  foot  from  the  junction  to  monitor  the 
reflected  and  transmitted  waveform. 

The  results  are  shown  in  Figure  21.  The  reflected  signal  added  to 
the  incident  signal  (voltage  reflection  coefficient  was  negative)  and 
the  transmitted  signal  increased  when  a  bend  was  introduced  into  the  line. 
The  reflection  is  small  (less  than  6  percent)  for  bends  up  to  90  degrees. 
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ANGLE  OF  BEND,  n  —  deg 


LA-7995-116 


FIGURE  21  REFLECTION  AND  TRANSMISSION  COEFFICIENT  OF  A  BEND  IN  A  SINGLE 
WIRE  LINE  OVER  GROUND  AS  A  FUNCTION  OF  BEND  ANGLE 


Similarly,  the  transmitted  signal  varied  only  a  few  percent  up  to 
9  =  90°.  At  angles  greater  than  90  degrees,  the  reflection  and  trans¬ 
mission  coefficients  increased  sharply.  Of  course,  for  0  >  90°,  the  wave 
on  the  incident  arm  can  couple  to  the  transmitted  arm  before  the  incident 
signal  reaches  the  bend. 

The  variation  of  results  with  height  was  small  and  within  the  data 


uncertainties. 


The  perturbance  due  to  the  bend  was  time  dependent.  The  values 
shown  in  Figure  21  are  peak  values.  Generally  the  effect  of  the  bend 
was  relaxed  within  a  few  ns  (or  equivalently,  a  few  line  heights) .  This 
result  is  consistent  with  King’s  analysis,  which  shows  that  the  effects 
on  £(w)  and  c(w)  are  strong  only  near  the  junction.  At  later  times — or 
longer  wavelengths — the  strong  effects  near  the  junction  are  smoothed, 
since  perturbances  that  extend  over  only  a  small  fraction  of  a  wavelength 
cannot  be  resolved.  The  results  of  King’s  analysis  are  compared  to  mea¬ 
sured  results  in  Figure  22.  Since  the  measured  and  calculated  results 
are  in  reasonable  agreement,  the  effect  of  the  line  height’s  comparability 
to  the  wavelength  must  be  small  for  the  rise  times  considered.  In  any 
case,  the  effect  is  small  for  9  s  90°. 

The  effect  of  radius  of  curvature  on  the  reflection  coefficient  was 
briefly  investigated.  A  line  of  0.01-inch  diameter  was  shaped  into  a 


BEND  ANGLE,  0  —  deg 
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FIGURE  22  COMPARISON  OF  CALCULATED  AND  MEASURED  REFLECTION  COEFFICIENTS 
OF  A  BEND  AS  A  FUNCTION  OF  BEND  ANGLE 
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‘lO-degree  bend  with  the  radius  of  curvature  varying  from  24  inches  to 
'^.01  inch.  The  results  are  shown  in  Figure  23.  The  reflection  coeffi¬ 
cient  decreases  monotonically  as  the  radius  of  curvature  is  increased, 

B.  Symmetrical  Junctions 

The  reflected  and  transmitted  signals  were  also  measured  at  symmetri¬ 
cal  junctions.  The  line  configuration  and  results  are  shown  in  Figure  24. 
The  line  height  used  was  two  feet,  or  six  times  the  equivalent  length  of 

the  rise  time  (L  =  ct  ) . 

eq  rise 

If  the  two  lines  forming  the  junction  are  considered  as  two  lines 

of  characteristic  impedance  in  parallel  and  are  fed  by  a  single  input 

line  of  Z  ,  then  the  reflection  coefficient  would  be  0.33  and  the  trans- 
o 

mission  coefficient  would  be  0.667.  These  values  are  indicated  in  Figure 
24.  The  results  show  that  this  model  accurately  reflects  the  interaction 


RADIUS  OF  CURVATURE.  R  —  in. 
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FIGURE  23  REFLECTION  COEFFICIENT  FOR  A  90°  BEND  AS  A  FUNCTION  OF 
RADIUS  OF  CURVATURE 
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JUNCTION  ANGLE,  0  —  deg 
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FIGURE  24  REFLECTION  AND  TRANSMISSION  COEFFICIENT  AT  A  SYMMETRICAL 
JUNCTION,  WHERE  h  =  2  FEET 

at  the  junction  when  the  junction  angle  is  large.  Ihis  result  is  rea¬ 
sonable  since  at  large  angles  the  junction  lines  are  relatively  indepen¬ 
dent.  However,  when  the  angle  is  small,  the  junction  lines  are  more 
tightly  coupled  to  each  other,  making  the  net  Impedance  closer  to  Z^. 
This  reduces  the  reflection  and  transmission  coefficients.  The  reflec¬ 
tion  coefficient  is  shown  for  early  times  <~i  ns)  and  for  late  times 
(~8.5  ns).  At  early  times  the  reflection  coefficient  is  much  lower 
than  would  be  expected  from  the  uncoupled  line  of  characteristic  imped¬ 
ance  Z  .  In  the  limit  as  9  0,  the  wire  diameter  roughly  doubles  and, 

o 

since  Z  is  proportional  to  in  2h/a,  the  value  of  Z  for  the  junction 
o  '  o 
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wire  almost  equals  the  value  of  Z  for  the  input  line,  result! 

o 

a  very  small  reflection.  This  is  shown  as  the  dashed  portion  ■ 
early  time  curve  from  0  =  O'*  to  0  =  22.5°. 

One  may  calculate  the  local  Z  of  the  junction  lines  cons 

o 

two  wires  over  ground.  The  Z  for  this  configuration  is" 

o 


r» 


o 


30  J0n 


i*} 

d 


k  ^  mV 


where  h  =  line  height,  d  =  line  diameter,  and  D  =  line  spacing 

diverging  junction  wires,  D  increases  as  the  distance  from  the 

increases.  As  D  -»  «>,  Z  reaches  1/2  that  of  the  input  line  an 

o 

flection  coefficient  reaches  0.33, 

Wo  have  computed  the  local  impedance  at  distances  corresp 

1  ns  and  8,5  ns  after  the  wave  reaches  the  junction.  The  refl 

coefficient  p  was  then  calculated  based  on  this  local  Z  as  a 

o 

the  input  line.  The  calculated  and  measured  results  are  shown 
following  tabulation: 


t  =  1  ns 

t  = 

8.5  ns 

0 

p  , 

0 

p  , 

0 

— 

calc 

meas 

calc 

meas 

22.5° 

0.192 

0.181 

0.285 

0.25 

45 

0,224 

0.234 

0.311 

0.283 

90 

0.264 

0.266 

0.328 

0.325 

135 

0.30 

0.334 

0.334 

0.342 

180 

0.333 

C.344 

0,333 

0.367 

ftie  calculation  does  not  include  the  effect  of  cumulativ> 
Lion  from  the  continuously  varying  impedance,  it  is  to  be  expe- 
the  values  would  not  be  in  perfect  agreement.  But  the  agreeme. 
reasonably  go.-.d  and  certainly  shows  the  correct  trends. 
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The  transmission  coefficient  may  be  calculated  from  the  reflection 
coefficient.  The  result  of  this  calculation  and  the  measured  values  are 


shown  in  the  following  tabulation: 


8_ 

calc 

neas 

22.5° 

0.596 

0.625 

45 

0.612 

0.625 

90 

0.632 

0.656 

135 

0.65 

0.656 

180 

0.666 

0.687 

Once  again,  the  agreement  is  reasonably  good,  and  the  trend  is  correct. 

Similar  measurements  were  carried  out  at  a  line  height  of  one  foct 
with  the  results  shown  in  Figure  25.  Calculations  as  outlined  above 
gave  similarly  close  agreement  with  the  measurements. 
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FIGURE  25  REFLECTION  COEFFICIENT  AT  A  SYMMETRICAL  JUNCTION.  WHERE 
h  =  1  FOOT 
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C,  Rise  Time  Characteristics 

The  waveform  of  the  incident  signal  is  shown  in  Figure  26.  The  rise 
time  is  310  ps.  For  0  of  22.5  degrees  and  45  degrees  the  rise  time  was 
measured  as  320  ps,  which  3  within  experimental  error.  As  6  was  increased 
beyond  45  degrees,  a  marked  degradation  in  rise  time  was  noted.  This  is 
illustrated  in  Figure  26  for  9  =  180°.  A  long,  slow  rise  is  observed. 

It  precedes  a  rapid  ncrease  of  signal  with  a  rise  time  approximately 
equal  to  the  input  rise  time.  The  difference  between  the  calculated 
time  it  would  take  a  signal  to  radiate  directly  from  the  feed  point  to 
the  observation  point  and  the  calculated  time  it  would  take  a  signal  to 
propagate  along  the  line  to  the  observation  point  equals  the  time  ahead 
of  the  fast  rise  portion  of  the  pulse  at  which  the  signal  is  first  ob¬ 
served.  We  conclude  that  the  slow  rising  "precursor"  signal  is  due  to 
coupling  directly  from  the  feed  point  radiation.  Even  when  the  feed 
point  is  infinitely  far  away,  the  electric  field  lines  associated  with 
the  transmission  line  mode  of  the  incident  wave  will  couple  to  the  junc¬ 
tion  arms  producing  a  current  flow  before  the  current  on  the  transmission 
line  flows  around  to  the  observation  point.  We  thus  expect  that  spur 
lines  will  see  such  a  "precursor"  even  when  the  signal  coupled  to  the 
main  line  is  vor,  far  from  the  ^pur. 
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FIGURE  26  CURRENT 
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SECTION  V 


COITPLING  TO  MULTI-WIRE  LINES 


Measurements  were  made  on  several  types  of  multi-wire  lines  for  both 
vertical  and  horizontal  polarization  and  for  a  range  of  azimuth  angles. 

A  sketch  of  the  three  wire  set-up  for  horizontal  polarization  is  shown 
in  Figure  27^  which  illustrates  some  of  the  nomenclature. 

The  phase  wires  are  labeled  1,  2,  and  3  while  the  neutral  is  labeled 
4.  The  height  of  the  phase  wires  above  ground  is  h  and  the  spacing  between 
wires  is  d.  The  electric  field  at  the  loads  was  54  V/m. 


The  azimuth  angle  is  cp.  Currents  were  monitored  at  points  A,  B, 
and  C  on  the  phase  wires.  This  point  was  about  4  feet  from  the  load  so 
that  there  was  sufficient  time  to  measure  the  currents  before  interaction 
with  the  load  occurred.  Current  at  point  D  (called  the  side-arm  current) 
gave  a  measure  of  the  differential  mode  current  that  would  flow  into  a 
load.  Measurements  at  point  E  monitored  the  current,  and  hence  potential, 
that  flowed  from  the  phase  to  the  neutral  wire.  As  described  later,  some 
measurements  were  performed  with  the  neutral  wire  periodically  grounded. 

To  gain  some  insight  into  the  complex  interactions  that  occur  in 
multi-wire  lines,  we  first  made  measurements  on  a  two-wire  line  where  the 
spacing  between  wires,  d,  was  2  feet  and  the  height  above  ground  was  also 
2  feet.  These  spacings  allowed  us  to  time  resolve  signals.  The  wave 
shapes  of  the  current  at  different  positions  on  the  line  are  shown  in 
Figure  28.  In  Figure  28(a)  the  current  is  shown  on  the  line  closest  to 
the  transmitted  signal.  It  rises  linearly  for  about  2  ns,  which  corre¬ 
sponds  to  the  time  before  the  ground-reflected  signal  arrives  at  the 
wire  (t  =  2h/c  sin  i|:) .  After  this  time  it  slowly  increases,  as  would  be 
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FIGURE  28  CURRENTS  COUPLED  TO  TWO-WIRE  LINES  OVER  GROUND 


expected  for  reflection  from  an  imperfect  ground.  The  measured  current 
at  2  ns  is  just  what  one  would  calculate  for  a  single  wire  over  ground. 


In  Figure  28(b)  the  current  is  shown  for  the  wire  that  is  farthest 
from  the  transmitter.  Its  waveform  is  similar  to  the  previous  one  except 
it  is  only  about  two-thirds  the  amplitude  at  a  given  time  after  the  ar¬ 
rival  of  the  wave.  This  reduction  in  amplitude  is  due  to  the  scattered 
fields  from  the  first  wire  inducing  a  current  that  opposes  the  current 
due  to  the  incident  field. 

Figure  28(c)  shows  the  current  at  the  short  circuit  load  at  a  point 
midway  between  the  wires.  The  current  reaches  a  peak  at  2  ns,  as  the 
current  from  the  first  wire  arrives  at  the  current  probe.  At  this  time, 
current  from  the  second  wire  arrives  and  opposes  the  current  to  a  reduced 
value.  The  current  does  not  go  to  zero  but  is  much  less  than  would  be 
collected  from  two  wires  in  free  space  as  analyzed  by  Harrison,®  Their 
analysis  shows  that  for  two  wires  in  free  space  the  differential  mode 
current  is 


E  d 

I  ,  =  -  cos  (40) 

d  Z 

o 

The  measured  value  in  Figure  28(c)  where  the  wires  are  over  ground, 
is  about  0.2  of  this  value.  We  believe  that  this  value  is  due  to  the 
ground  reflection  coefficient’s  being  about  -0.8  instead  of  -1  as  it 
would  be  for  perfect  ground.  Our  reasoning  is  explained  as  follows. 

The  differential  current  flowing  in  the  short-circuit  load  of  a 
two-wire  line  in  free  space  is  as  given  in  Eq.  (40).  Placing  these 
wires  over  ground  is  approximately  equivalent  to  exciting  these  wire 
from  an  angle,  -  \jr,  and  with  a  field  strength  of  RE  ,  where  R  is  the 


ground  reflection  coefficient.  Therefore,  the  differential  current  due 
to  presence  of  ground  is 


t 

dg 


RE  d  cos  (-t^) 
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RE  d  cos  ti 
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Z 
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(41) 


Since  R  is  a  negative  number  less  than  one,  the  net  differential  current 
after  the  ground  reflection  will  be 


E  d  cos 

I  =  (l  -  -  (42) 

d  net  '  '  Z 

o 

This  equation  is  only  approximately  true  in  this  case  since  the  presence 
of  ground  will  change  Z^  somewhat.  Nevertheless,  the  predominant  effect 
of  ground  will  be  to  reduce  the  differential  current. 

Measurements  on  four-wire  lines  consisting  of  three  phase  wires  and 
a  neutral  were  made  with  the  results  shown  in  Figure  29.  The  wires  were 
terminated  in  400  ohm  loads.  Once  again  the  conductor  closest  to  the 
transmitter  had  the  largest  current  flowing  in  it.  Figure  29  shows  the 
current  in  the  side  arm,  which  would  produce  differential  voltages.  The 
peaks  values  are  about  10  percent  of  the  current  flowing  on  the  lines 
while  the  late  time  currents  are  only  a  few  percent  of  the  line  currents. 

For  this  20:1  model  of  a  power  line,  the  full  scale  currents  are 
shown  in  parentheses  for  an  incident  field  of  54  V/m.  If  the  incident 
field  is  higher  the  current  will  also  be  higher  by  the  ratio  of  the  fields. 

The  phase  to  neutral  currents  are  shown  in  Figure  30.  These  currents, 
even  at  late  times  are  at  least  10  percent  of  the  line  currents. 

The  same  measurements  were  made  with  the  neutral  grounded  every  2 
feet.  The  currents  are  shown  in  Figure  31.  The  line  current  has  not 
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HORIZONTAL  POLARIZATION 
FOUR-WIRE  LINE 
400  n  LOADS 
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FIGURE  30  PHASE-NEUTRAL  CURRENTS  ON  A  FOUR-WIRE  LINE 
WITH  400  n  LOADS  AND  HORIZONTAL  POLARIZATION 

changed  much  from  the  line  current  in  the  ungrounded  case.  However,  the 
late  time  and  the  phase  to  neutral  current  has  increased  significantly 
over  the  ungrounded  case.  The  direction  of  current  is  such  that  current 
is  being  driven  from  the  neutral  into  the  phase  wires. 

The  four-wire  system  with  multiply  grounded  neutral  was  also  excited 
with  vertical  polarization.  The  results  are  shown  in  Figure  32.  The 
line  current  is  a  maximum  at  0  degrees  and  dips  to  a  minimum  at  90  degrees 
where  the  wires  are  cross  polarized  to  the  wave.  The  side  arm  currents 
at  0  degrees  are  reduced  to  about  10  percent  of  the  line  current.  The 
phase  to  neutral  currents,  however,  are  quite  large,  being  equal  to  the 
line  currents  at  0  degrees.  They  are  larger  by  an  order  of  magnitude  at 
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FIGURE  32  CURRENTS  COUPLPO  TO  A  FOUR-WIRE  LINE  WITH  400  fi  LOADS, 
MULTIPLY  GROUNDED  NEUTRAL,  AND  VERTICAL  POLARIZATION 


90  degrees  because  the  vertical  grounding  elements  pick  up  current  at 
90  degrees. 

Thus  the  differential  currents  in  ungrounded  systems  are  typically 
an  order  of  magnitude  below  the  line  currents.  However,  in  grounded 
neutral  systems,  large  currents,  comparable  to  the  line  cuxrents,  can 
flow  frcm  the  neutral  to  the  phase  wires. 


SECTION  VI 


CONCLUSIONS  AND  RECONDtENDATIONS 

As  a  result  of  the  work  described  in  this  report  and  in  earlier 
documents  prepared  under  this  task,^  it  is  believed  that  the  coupling 
mechanisms  by  which  the  EMP  interacts  with  open  transmission  lines  such 
as  power  distribution  lines  is  well  understood  and  is  adequately  described 
for  engineering  purposes  by  the  low  frequency  analysis  used  heretofore. 

It  is  apparent  from  the  theoretical  analysis  and  experimental  verifica¬ 
tion  of  this  analysis,  that  the  strongest  common-mode  coupling  to  the 
power  distribution  lines  occurs  for  a  vertically  polarized  wave  incident 
almost  along  the  axis  of  the  transmission  line.  For  this  rather  special 
case,  the  short-circuit  current  induced  in  a  straight  line  1  kilometer 
long  is  of  the  order  of  1  ampere  for  an  incident  field  strength  of  I  volt 
per  neter. 

In  addition  to  verifying  the  coupling  theory,  the  work  reported  here 
demonstrated  the  effects  of  several  common  perturbations  to  the  ideal 
straight  transmission  line.  These  include  bends  and  junctions,  as  well 
as  periodic  grounds,  and  differential  mode  coupling.  The  effects  of  a 
bend  in  the  transmission  line  were  found  to  be  measureable  and  interest¬ 
ing,  but  of  little  significance  in  terras  of  propagation  of  the  E3i!P-induced 
common  mode  current  except  when  the  bend  angle  approaches  180  degrees. 
Likewise,  the  effects  of  junctions  of  two  lines  with  a  third  were,  for 
most  cases,  found  to  be  adequately  described  by  conventional  transmission- 
line  theory,  although  deviations  from  ideal  transmission-line  behavior 
were  noted  when  the  two  load  lines  were  quite  close  together.  Differen¬ 
tial  mode  currents,  on  the  other  hand,  could  be  significant  since  they 
can  be  as  much  as  1  percent  of  the  maxiraum  common-node  current.  The 
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mari^uin  differential  current  occurs  for  broadside  incidence,  however, 
v/hereas  the  uaxi-riir  conunon-mode  current  occurs  for  end-on  incidence. 

One  of  xlie  results  of  this  study  was  to  have  been  recomendations 
for  full-scale  tests  to  obtain  additional  information  on  coupling  to 
power  lines  or  propagation  on  power  lines.  It  is  felt,  however,  that 
the  results  of  the  scale- model  experiments  and  analysis  have  adequately 
established  the  coupling  and  propagation  characteristics  of  typical 
transnissio.-line  configurations.  Furthermore,  the  return  from  full- 
scale  ‘.ests  is  so  low  (in  relation  to  cost)  that  such  testing  should  be 
reserved  to  testing  related  to  weapon  systems,  in  which  the  interaction 
of  the  existing  power  system  with  the  weapon  system  and  the  local  terrain 
is  to  be  determired.  There  are  few  aspects  of  the  physics  of  cottpling 
to,  and  propagation  along,  power  lines  that  catinot  be  determined  more 
efficiently  under  controlled  laboratory  conditions  than  under  field-test 
conditions.  Parameters  such  as  soil  conductivity  and  stratification  or 
local  terrain  in  the  vicinity  of  a  particular  system  are  not  easily 
assessed  in  the  general  scale-model  study;  however,  these  peculiarities 
of  a  specific  system  or  site  could  n'kt  be  evaluated  from  an  arbitrary  full- 
scale  test,  either.  Only  a  foil-scale  test  on  the  site  of  interest  would 
provide  the  information  necessary  to  assess  these  site- peculiar  aspects 
of  the  power  system.  Therefore,  full-scale  tests  for  the  purpose  of 
evaluating  coupling  and  propagation  charactezistlCF  are  not  recommended, 
although  continued  testing  of  the  vc-apon  system  res|ionse  co  the  coupled 
signal  will  be  necessary. 
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Appendix  A 

COMPARISON  OF  TRANSMISSION-LINE  APIUOXIMATION 
WITH  THE  EXACT  SOUJTION 


The  transmission-line  analysis  of  the  current  induced  in  the  poser 
line  is  based  on  the  assumption  that  the  characteristic  impedance  of  the 
line  remains  constant  at  the  value  defined  by  the  electrostatic  capaci¬ 
tance  pex*  unit  length  and  the  magnetostatic  inductance  per  unit  length. 
When  the  line  height  is  comparable  to  or  greater  than  a  wavelength,  the 
characteristic  impedance  so  defined  is  no  longer  valid.  As  has  been 
shown  in  this  report,  however,  the  deviation  of  predicted  values  from 
the  measured  values  is  so  small  that,  within  the  accuracy  of  the  measure¬ 
ment,  it  is  not  detectable.  A  more  rigorous  analysis  of  the  current  in¬ 
duced  in  a  perfectly  conducting  cylinder  over  a  perfectly  conducting 
groimd  plane  has  recently  been  made  by  Flammer  and  Singhaus.*^  The 
purpose  of  this  appendix  is  to  compare  the  results  of  the  transmission- 
line  analysis  with  the  exact  solutions  obtained  in  Ref.  10. 

The  comparison  will  be  made  for  a  conductor  height-to-radius  ratio 
h/a  =  100,  which  coriesponds  to  a  transmi£.>ion  line  whose  effective 
radius  is  10  cm  at  a  height  of  10  m.  The  characteristic  impedance  of 
this  line  is  317  ohms.  The  transmission-line  solutions  for  the  step- 
function  response  of  the  current  are,  for  vertical  polarization. 
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The  characteristic  impedance  has  been  expressed  as 
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to  accommodate  the  current  normalization  used  by  Flammer  and  Singhaus. 
Klammer  and  Singhaus  also  use  a  dimensionless  time-space  parameter 

t 

ct  -  z  COSY 

u  =  - 

a  sinY 


which,  when  the  spatial  variation  is  suppressed,  becomes 


ct 

_  - -  (for  vertical  polarization,  cp  =  0) 

V  a  sini^ 


or 


(for  horizontal  polarization,  co  =  90°) 


(The  angle  7  used  by  Flammer  and  Singhaus  is  the  polar  angle  measured 
irom  the  axis  of  the  wire;  see  Figure  1  for  the  elevation  and  azimuth 
convention  used  in  this  report.)  In  the  examples  that  follow,  time  will 
be  xpressed  in  graduations  of  (h  sin\^)/c  or  h/c  instead  of  the  dimension- 
iGos  parameter  u. 

The  transmission-line  solution  and  the  exact  solution  for  the  cur¬ 
rent  induced  in  an  infinitely  long  conductor  by  a  step  function  of  inci¬ 
dent  field  are  shown  in  Figures  A-1  and  A-2.  Figure  A-1  shows  the  re- 
onsf’  for  a  vertically  polarized  incident  field  for  an  azimuth  angle  of 
idence  of  fp  =  0  and  arbitrary  elevation  angle  of  incidence.  Figure  A-2 
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FIGURE  A-1  COMPARISON  OF  THE  TRANSMISSION-LINE  APPROXIMATION  AND  THE 
EXACT  SOLUTION  FOR  THE  CURRENT  IN  A  WIRE  OVER  A  PERFECTLY 
CONDUCTING  GROUND  PLANE  —  STEP-FUNCTION  INCIDENT  FIELD. 
VERTICALLY  POLARIZED 


shows  the  response  for  a  horizontally  polarized  incident  field  at  an 
azimuth  angle  of  90°  (broadside)  for  selected  elevation  angles  of  inci¬ 
dence.  It  is  apparent  from  these  illustrations  that  the  transmission¬ 
line  approximation  is  a  good  approximation  to  the  current  Induced  in  a 
wire  over  a  ground  plane  even  for  the  zero-rise-time  step  function.  For 
realizable  finite  rise-time  pulses^  the  difference  between  the  exact 
solution  and  the  transmission-line  approximation  should  be  so  small  that 


! 
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it  would  be  difficult  to  detect  experimentally  (as  has  been  observed). 

The  exact  solution  does  provide  interesting  insight  into  the  behavior  of 
the  scattered  waves  before  the  structure  "settles  down"  to  a  transmission 

I 

line.  In  Figure  A-2,  for  example,  one  can  see  the  discontinuity  that 
occurs  as  the  ground-reflected  wave  arrives  at  (2h  sin')()/c,  a  second 
discontinuity  when  the  scattered  direct  wave  reflected  from  the  ground 
returns  at  2h/c,  and  a  third  discontinuity  when  the  ground-reflected  wave 
scattered  from  the  wire  and  reflected  back  from  the  ground  arrives  at 
(1  +  sin-;)  2h/c.  Subsequent  multiply-scattered  waves  have  a  very  weak 
influence,  but  they  eventually  bring  the  exact  Zate-time  I’esponse  into 
coincidence  with  the  transmission-line  solution. 

A  comparison  of  the  impulse  responses  for  a  vertically  polarized  wave 
incident  at  ro  =  0  is  shown  in  Figure  .A-3.  The  zero-width  impulse  induces 
a  very  large  spike  when  the  direct  wave  arrives,  but  the  current  quickly 
decays  to  a  value  comparable  to  that  predicted  by  the  transmission-line 
approximation.  Since  the  zero-width  impulse  is  not  realizable,  the  dif¬ 
ference  between  the  exact  and  approximate  solutions  for  practical  short 
pulses  of  finite  width  and  height  and  of  non-zero  rise  and  fall  times  is 
probably  not  very  significant.  The  exact  solution  shows  a  small  follow- 


on  current  not  predicted  by  the  transmission- line  theory. 

A  lurther  comparison  is  shown  in  Figure  A-4,  where  the  response  to 
a  vertically  polarized  incident  exponential  pulse  is  shown.  This  wave- 
lonti  was  obtai.ied  by  the  convolution  of  an  exponential  waveform 


E(t)  =  Eq  e 


where 


20  h  siniii 


with  the  step  function  responses  of  Figure  A-1.  Because  the  time  constant 
t^  of  the  exponential  waveform  was  chosen  for  ease  of  computing,  it  is  a 
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FIGURE  A-3  COMPARISON  OF  THE  TRANSMISSION-LINE  APPROXIMATION  AND  THE 
EXACT  SOLUTION  FOR  THE  CURRENT  IN  A  WIRE  OVER  A  PERFECTLY 
CONDUCTING  GROUND  PLANE  —  IMPUlSE  INCIDENT  FIELD, 
VERTICALLY  POLARIZED 
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FIGURE  A-4  COMPARISON.  OF  THE  TRANSMISSION-LINE  APPROXIMATION  AND  THE 
EXACT  SOLUTION  FOR  A  VERTICALLY  POLARIZED  INCIDENT 
EXPONENTIAL  PULSE 

function  of  the  angle  of  incidence  for  this  example.  It  is  ten  times  the 
"clear  time"  between  the  beginning  of  the  pulse  and  the  arrival  of  the 
ground-reflected  wave,  so  that  for  a  10-m-high  line  with  the  wave  inci¬ 
dent  from  directly  overhead,  t^  =  610  ns,  whereas  for  incidence  at  a  30° 
elevation  angle,  t^  =  305  ns.  The  exponential  pulse  has  zero  rise  time. 

A  comparison  of  the  frequency  spectra  of  the  exact  solution  and  the 
transmission-line  approximation  is  shown  for  a  horizontally  polarized 
wave  incident  at  an  elevation  angle  of  45°  in  Figure  A-5.  This  compari¬ 
son  of  the  current  spectra  is  consistent  with  the  time-domain  comparisons 
in  that  the  low-frequency  spectra,  corresponding  to  the  late-time  re- 
spon.ce^  are  similar,  while  the  high-frequency  spectra,  corresponding  to 
the  early-time  response,  differ  somewhat.  Both  spectra  are  characterized 
by  similar  maxima  and  zeros. 
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Appendix  B 


THEORY  OF  ELECTROMAGNETIC  SCALING 


The  theory  of  models  is  widely  used  in  many  branches  of  engineering, 
especially  where  the  theoretical  difficulties  make  a  problem  intractable 
to  analytical  solution.  Before  discussing  the  modeling  relations  for 
Maxwell's  equation,  a  discussion  of  modeling  in  the  context  of  the  fa¬ 
miliar  RC  circuit  will  illustrate  some  of  the  features  of  modeling.  The 
differential  equation  for  this  circuit  is 

t 

V  -  M  +  ^  j  idt  (B-l) 

o 


Di f Ccrcntiating  with  respect  to  time,  and  rearranging,  we  obtain 


^  a-i- 
5t  ‘  RC 


0 


(B-2> 


Now  suppose  that  we  could  not  solve  this  equation,  but  wanted  to  make 
measurements  on  this  circuit  to  find  the  variation  of  current  with  re¬ 
spect  to  time  and  the  values  of  R  and  C.  This  is  analogous  to  our  being 
unable,  without  great  difficulty,  of  solving  for  the  current  as  a  func¬ 
tion  ol  time  in  our  transmission  line  geometry. 

I  1  we  had  the  proper  R  and  C  and  V,  we  could  simply  measure  the 
response.  However,  if  it  is  impractical  to  make  measurements  in  a  labora¬ 
tory  (say,  because  C  is  so  large),  we  could  still  study  the  response  by 
using  a  model  and  the  full-scale  system  can  be  derived  as  follows.  As¬ 
sume  the  following  relations  between  the  parameters  in  the  two  systems; 
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i  =  ki' 
t  =  mt' 

R  =  nR' 

C  =  pC'  (B-3) 


The  primed  quantities  refer  to  the  model  system  and  the  unprimed 
quantities  to  the  full-scale  system.  The  dimensionless  numbers  m,  n, 
and  p  are  the  scale  factors  for  the  current,  time,  resistance,  and  ca¬ 
pacitance,  respectively. 


For  the  model  system  the  differential  equation  is 


di' 

bt' 


R  C 


0 


(B-4) 


If  we  wish  to  make  the  current-time  response  similar  for  ooth  sys¬ 
tems,  then  we  should  be  able  to  apply  the  scaling  relations  of  Eq.  (B-3) 
to  Eq.  (B-4)  and  thus  show  under  what  conditions  Eq.  (B-2)  and  Eq.  (B-4) 
are  equal. 

For  i^  we  substitute  i/k.  Similarly,  with  r'  and  C'  we  substitute 
R/n  and  C/p.  For  biVbf  we  use 


_5 ^  ^ 

bt'  bt  bt'  bt 


(B-5) 


Equation  (B-4),  with  these  substitutions,  becomes 

^  ilE 

k  bt  k  RC  " 


This  equation  will  be  identical  to  Eq.  (B-2)  if 


np 

k 


1 


(B-6) 
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These  define  the  scaling  relations  for  similitude. 

The  first  condition  says  that  the  current  is  scaled  by  the  same 
factor  as  is  the  time,  if  the  time  is  compressed  by  a  factor  of  100, 
the  model  system  current  will  be  decreased  by  the  same  factor.  Full- 
scale  currents  can  be  found  from  model -measured  currents  by  multiplying 
by  m. 

The  second  condition  says  that  the  product  of  the  resistance  and 
capacitance  scale  factors  must  equal  the  current  scale  factor.  Since  the 
current  and  time  factors  are  the  same  from  the  first  condition,  this  is 
the  same  as  saying  that  the  resistance  and  capacitance  is  changed  in  the 
same  proportion  as  the  time.  If  time  is  compressed  by  a  factor  of  100, 
the  product  of  resistance  and  capacitance  must  be  decreased  by  the  same 
factor. 

This,  of  course,  is  just  restating  the  result  known  from  the  solution 
oI  Kq.  (B-2)  that  the  current  depends  on  the  ratio  t/RC  and  not  t  or  RC 
separately.  Note  that  this  conclusion  may  be  derived  from  an  analysis 
of  the  scaling  relations  without  actually  solving  the  equation. 

This  procedure  can  be  applied  to  Maxwell's  equations  to  find  the 
scaling  relations  for  electromagnetic  problems.  The  conditions  for 
similitude  will  lead  to  relations  between  the  parameters  in  Maxwell's 
equations.  When  these  conditions  hold,  the  relation  between  the  full 
scale  fields  and  parameters  will  be  similar  to  those  for  the  model  system. 
The  scaling  relations  will  enable  one  to  apply  results  from  the  model 
.system  to  the  full  scale  system. 

The  most  complete  discussion  of  electromagnetic  scaling  i?  given  by 
Sinclair.’^  We  shall  draw  on  his  discussion  for  most  of  what  follows. 
Consider  a  full  scale  system  that  is  described  by  spatial  coordinates 
X,  y,  time  coordinate,  t,  and  el<ctric  and  magnetic  field  vectors 
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E  vx,  y,  z,  t),  and  H  (x,  y,  z,  t).  The  corresponding  quantities  in  the 
model  system  are  primed.  Sinclair  shoAs  that  if  there  are  scale  factors 
that  relate  the  full  scale  to  the  model  system  quantities  such  that 
X  =  px' 

y  =  py' 
z  =  pz* 
t  =  'rt' 


E  (X,  y,  z. 

t)  =  oE' 

(x'. 

/  /  .  /  v 

y  ,  ^  ^ 

H  (X,  y,  z. 

t)  =  SH' 

(x'. 

y  ^  2  ,  t  ) 

a  set  of  relations  between  the  scale  factors  (p,  Yj  O'-  0)  may  be  found 
that  ensure  that  the  model  is  an  accurate  simulation  of  the  full-scale 
.system.  This  is  done  by  equating  Maxwell's  equations  in  the  full-scale 
system  to  a  set  of  equations  derived  from  Maxwell’s  equation  in  the  model 
system  that  are  transformed  to  full-scale  quantities  by  the  appropriate 
scale  factors. 

Thus  Ma.xwell's  equations  in  the  full-scale  system  are 

bE 

curl  H  =  cE  -r  €  — 

Ot 

curl  E  =  -►.  —  (B-8) 

ot 


In  the  model  system  they  are 
curl^  = 

curl^  e'  = 

To  illustrate  the  procedure, 

curl'  ,  H'  = 

X 


Ml 

bt' 


consider  the  term 

bH',  bH', 

_ 5 _ ^ 

by'  bz' 


(B-9) 


(B-10) 
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Since  H  =  8H'  and 


6  b  by  ^ 

5y'  by  by'  ^  by 

curl',  H'  =  ^  curl  H  (B-11) 

X  e 


In  this  way  each  of  the  model  system  terms  are  related  to  the  full- 
scale  system. 


A  set  of  three  equations  that  must  be  satisfied  ensures  accurate 
simulation.  These  are 


P  3 

- ji 

Y  O' 


(B-12) 


Sinclair's  article  continues  with  the  general  scaling  properties,  but  for 
our  casc.s  we  can  deduce  several  important  conclusions  from  the  above  set 
of  equations.  If  permittivity  and  permeability  arc  to  be  the  same  in 
both  the  model  and  the  full-scale  system,  then  cr/g  =  i,  p/Y  =;  1,  and 
n'  =  pc.  The  first  condition  says  that  the  ratio  of  E/H  in  the  two  sys¬ 
tems  is  the  same.  This  is  satisfied  if  the  e  and  u  in  both  systems  are 
the  same.  The  second  condition  requires  that  the  length  scale  and  the 
time  scale  must  vary  by  the  same  factor.  If  the  model  dimensions  are 
1 '100  of  full  scale,  then  the  time  must  also  be  compressed  by  a  factor 
of  100. 

The  third  condition  says  that  the  conductivity  in  the  model  system 
must  be  changed  directly  with  the  length  scale  factor.  The  mode’  con¬ 
ductivity  must  be  100  times  the  full-scale  conductivity  if  the  model  di¬ 
mensions  are  1/100  of  the  full  scale  dimension. 
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If  the  above  relations  are  satisfied,  then  many  other  relations  can 
be  derived,  such  as  Poynting  vector,  power,  gain,  echo  area,  current 
density,  etc.  The  quantity  most  important  for  us  is  that  for  current: 


1  =  Sp  I 


(B-13) 


The  full-scale  current  is  equal  to  the  model  current  multiplied  by 
the  length  scale  factor  and  the  ratio  of  the  full-scale  magnetic  field 
to  the  model  magnetic  field.  Since  o  =  8  we  could  as  well  use  the  ratio 
of  electric  fields. 


I  =  o-pl 


(B-14) 


If  one  multiplies  and  divides  by  E,  one  obtains 


1  =E  pp 


(B-15) 


In  our  measurements  we  will  determine  l' /e' ,  so  that  full-scale 
currents  may  be  found  by  multiplying  by  Ep. 

The  restrictions  on  such  modeling  require  that  the  media  be  linear. 
The  media  may  vary  from  point  to  point  in  space  but  mu.st  be  independent 
of  tine. 


A  practical  restriction  that  is  usually  unimportant  is  that  if  the 
full-scale  system  has  conductors  (like  antenna  elements)  as  well  as  a 
conducting  ground,  the  conductivity  of  the  conductors  should  be  modeled 
for  perfect  simulation.  This  is  dif«'icoi'  do  if  the  full-scale  con¬ 
ductors  are  copper.  However,  in  cases  :re  the  losses  on  the  conductors 
are  negligible,  the  effect  of  not  modeling  the  conductor  conductivity  is 
unirojwrtant.  For  example,  with  a  rod-over-ground  transmission  line  it 
is  usually  unimp<'rtant  to  scale  conductivity  of  the  rod,  since  most  of 
the  loss  is  in  the  ground.  Only  the  ground  conductivity  must  then  be 
s-ralcd  for  proper  simulation. 
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Appendix  C 


ESTIMATED  ACCURACY  OF  THE  MODEL  DATA 


The  scale-model  experiments  were  conducted  with  a  pulsed  dipole  an¬ 
tenna  illuminating  the  models  and  measurements  of  the  response  of  the 
model  were  made  with  a  current  probe  and  recorded  with  an  oscilloscope 
and  X-Y  recorder  as  illustrated  in  Figure  C-1.  The  components  of  this 
system  that  are  important  in  evaluating  the  accuracy  of  the  measurements 
are  the  current  probe^  the  oscilloscope^  and  the  transmitting  antenna. 

The  current  on  the  transmission  line  was  measured  with  a  Tektronix 
CT-1.  This  current  probe  has  an  accuracy  of  ±5  percent.  The  currents 
were  detected  using  a  Hewlett-Packard  sampling-scope  system  with  an 
accuracy  of  ±3  percent. 

The  largest  error  is  that  associated  with  the  peak  radiated  field. 
There  are  frequency-dependent  losses  in  the  cable  and  balun  that  feed  the 
dipole  radiator.  A  measurement  of  the  efficiency  of  conversion  of  dc 
supply  voltage  to  current  at  the  antenna  terminals  gave  a  value  of  about 
70  percen’c  of  the  theoretical  value  when  no  losses  were  assumed. 

The  dipole  had  tapered  tips  at  the  feedpoint  for  4  inches  of  length 
and  then  was  a  uniform  cylinder  of  1  inch  diameter.  This  structure  was 
treated  as  a  uni foxm  bicone  of  diameter  equal  to  1  inch  and  length  as¬ 
sociated  with  the  ground  clear  time  (2h/c  sin^)  to  evaluate  the  peak 
field.  It  is  estimated  that  the  peak  field  per  dc  charge  volt  could  be 
in  error  by  ±15  percent. 

The  1/e  decay  time  of  the  current  at  the  antenna  feedpoint  is  35  ns. 


This  represents  a  decrease  in  field  of  about  4  nercent  in  2h/c  sint  for 


SALTED  SOIL  (O  0.2-0.3  mho/m) 

la-7995-133 


FIGURE  C-1  SCALE-MODEL  EXPERIMENT  TO  STUDY  COUPLING  TO  TRANSMISSION  LINE 

the  2-foot-high  line,  and  was  accounted  for  in  Ref.  2.  At  much  later 
times,  the  waveform  will  have  larger  deviations  from  a  step-function  input 
that  will  also  include  a  decrease  in  radiated  field  due  Lo  increasing 
impedance  of  the  cylindrical  dipole  antenna. 

Because  the  radiating  source  was  a  finite  distance  from  the'  power¬ 
line  model,  different  parts  of  the  line  were  illuminated  with  different 
amplitudes  and  polarizations.  For  co  =  90®  illumination,  the  errors  were 
quite  small  (a  few  percent)  for  the  peak  values,  increasing  slowly  with 
time.  At  cp  close  to  0°  the  errors  were  larger,  although  even  there,  the 
peak  field  error  was  estimated  to  be  no  more  than  about  10  percent. 

Therefore  if  the  current  probe  and  the  oscilloscope  errors  are  as¬ 
sumed  to  be  random  so  that  they  combine  in  a  root-mean-square  fashion, 
the  random  error  in  the  measurements  will  be  about  6  percent.  In  addi¬ 
tion,  there  are  accountable  systematic  deviations  due  to  the  droop  in 
the  radiated  "step  function"  and  the  nonuniformity  of  the  field.  Finally, 
there  is  an  uncertainty,  which  is  probably  systematic,  of  about  15  per¬ 
cent  in  absolute  magnitude  of  the  incident  field  strength. 
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